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Abstract We introduce a new class of distances between nonnegative Radon measures in 
M. d . They are modeled on the dynamical characterization of the Kantorovich-Rubinstein- 
Wasserstein distances proposed by BENAMOU-BRENIER |7| and provide a wide family in- 
terpolating between the Wasserstein and the homogeneous W y lp -Sobolev distances. 
From the point of view of optimal transport theory, these distances minimize a dynamical 
cost to move a given initial distribution of mass to a final configuration. An important dif- 
ference with the classical setting in mass transport theory is that the cost not only depends 
on the velocity of the moving particles but also on the densities of the intermediate configu- 
rations with respect to a given reference measure y. 

We study the topological and geometric properties of these new distances, comparing them 
with the notion of weak convergence of measures and the well established Kantorovich- 
Rubinstein- Wasserstein theory. An example of possible applications to the geometric theory 
of gradient flows is also given. 

Keywords Optimal transport ■ Kantorovich-Rubinstein- Wasserstein distance ■ Continuity 
equation ■ Gradient flows 



1 Introduction 

Starting from the contributions by Y. Brenier, R. McCann, W. Gangbo, L.C. Evans, 
F. Otto, C. Villani (5, 18,25, 17,27], the theory of Optimal Transportation has received 
a lot of attention and many deep applications to various mathematical fields, such as PDE's, 
Calculus of Variations, functional and geometric inequalities, geometry of metric-measure 
spaces, have been found (we refer here to the monographs |28, 16,30,3,31]). Among all pos- 
sible transportation costs, those inducing the so-called ZAKANTOROVICH-RUBINSTEIN- 
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WASSERSTEIN distances W p (Ho,fl\),p&(l, +°°), between two probability measures fl,V6 



Wpi^fr) :=inf\ / b-x^di: : rer(/io,Mi)^ (i.i) 

play a distinguished role. Here F(Mo, Mi ) ls me set °f al l couplings between Mo and Mi : they 
are probability measures E on R d x R d whose first and second marginals are respectively 
Mo and Mi , i.e. E(B x R d ) = Mo(B) and E(R d x S) = Mi (B) for all Borel sets B e 

It was one of the most surprising achievement of [24 ,25 .19,26] that many evolution 
partial differential equations of the form 

d,p+V.(p|€r 2 S)=0, S=- V (^) inKM0,+~), (1-2) 

can be, at least formally, interpreted as gradient flows of suitable integral functionals & 
with respect to W p (see also the general approach developed in 1 30 3 ,31 1). In d 1 -2b / 8p 
is the Euler first variation of & , q := p/(p — 1) is the Holder's conjugate exponent of p, 
and 1 1— > p, (a time dependent solution of i ll. 2b ) can be interpreted as a flow of probability 
measures ji, = p, Jz? with density p, with respect to the Lebesgue measure Jt? d in R d . 

Besides showing deep relations with entropy estimates and functional inequalities [27|, 
this point of view provides a powerful variational method to prove existence of solutions to 
l ll.2l ). by the so-called Minimizing movement scheme 1 1 9 , 1 3 , 3 1 : given a time step X > and 
an initial datum Mo = Po-^ rf . the solution ji, = p,«Sf d at time t « m can be approximated by 
the discrete solution \i" obtained by a recursive minimization of the functional 

M^-^T^/(M,^)+^(M), * = 0,1,... (1.3) 

The link between the Wasserstein distance and equations exhibiting the characteristic struc- 
ture of ( 11.2b (in particular the presence of the diffusion coefficient p, the fact that £ is a 
gradient vector field, and the presence of the ^-duality map £j i— > |{j | 9 ~ 2 |), is well explained 
by the dynamic characterization of W p introduced by BENAMOU-BRENIER [7]: it relies in 
the minimization of the "action" integral functional 



W/(Mo,Mi)=inf{ j f ft(*)|v,(x)|"dx(te : 

F I- JO JR d 

<3,p,+V-(p,v,)=0inR d x(0,l), mo = P|, =0 ^, Mi=P|, =1 ^}- 



(1.4) 



Towards more general cost functionals. If one is interested to study the more general class 
of diffusion equations 

d t p+V-(h(p)\^\ q - 2 ^=0, |=-v(^) mR d x(0,+~), (1.5) 

obtained from d 1 -2b replacing the mobility coefficient p by an increasing nonlinear function 
h(p), h : [0, +°°) — ► [0, +°°) whose typical examples are the functions h(p) = p a , a > 0, it 
is then natural to investigate the properties of the "distance" 

W/(Mo,Mi)=inf{ C [ h(p t (x))\y,(x)\P±xdt: 

<9,p,+V.(/;(p,)v,)=0inR<'x(0,l), Mo = P|, =() ^, Mi = P| /=1 ^}- 
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In the limiting case a = 0, h(p) = 1, one can easily recognize that l ll.6l > provides an equiv- 
alent description of the homogeneous (dual) W- hp (R d ) Sobolev (pseudo)-distance 

\\to-frh-hP( M *)--™p{[jd(w-fr):£eCl(R d ), £jD£|«dx<l}. (1.7) 



Thus the distances defined by dl.6l > for < a < 1 (we shall see that this is the natural 
range for the parameter a) can be considered as a natural "interpolating" family between 
the Wasserstein and the (dual) Sobolev ones. 

Notice that if one wants to keep the usual transport interpretation given by a "dynamic 
cost" to be minimized along the solution of the continuity equation, one can simply introduce 
the velocity vector field \, := p,h(p,)\, and minimize the cost 

_p/(p) |f,|'cUdf, where f(p) := (^j)^ ( L 8) 

Therefore, in this model the usual p-energy p ( |v, \ p dx of the moving masses p, with 
velocity \ t results locally modified by a factor f(p t ) depending on the local density of the 
mass occupied at the time t . Different non-local models have been considered in |U|4j. 

In the present paper we try to present a systematic study of these families of intermediate 
distances, in view of possible applications, e.g., to the study of evolution equations like l ll.5l >. 
the Minimizing movement approach tl.3b . and functional inequalities. 

Examples: PDE's as gradient flows. Let us show a few examples evolution equations which 
can be formally interpreted as gradient flows of suitable integral functionals in this setting: 
the scalar conservation law 

d t p — V- (p a W) =0 corresponds to the linear functional &(p) ■= / V(x)pdx, 
for some smooth potential V : R d — > R and p = 2. Choosing for in > 



p = 2, &(p)=c a . m lp m+1 - a dx, 



m 



c, 



a.m 



(m + 1 — a)(m — a) ' 
one gets the porous media/fast diffusion equation 



d,p - -^—V ■ (p a Vp"'- a ) = d,p - Ap m = 0, (1.9) 
m — a 

and in particular the heat equation for the entropy functional ( 2 -a)\i-a) / P 2 " Choosing 

r m+iq-i-a m{q-\) q 

& P) = c a ,m.q / P d>c, c amq := - — — ' -, 

H J H (m + 2q-3-a)(m + q-2-a) 

one obtains the doubly nonlinear equation 

d,p - mV ■ (p"'- 1 \Vp\"- 2 Vp) = (1.10) 

and in particular the evolution equation for the g-Laplacian when m = 1. The Dirichlet 
integral for p = 2 

JF(p) = -J \Vp\ 2 dx yields d,p + V ■ (p a VAp) = 0, (1.11) 
a thin-film like equation. 
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The measure-theoretic point of view: Wasserstein distance. We present now the main points 
of our approach (see also, in a different context, 1101 ). First of all, even if the language of 
densities and vector fields (as p and \,\ in 1 11. 4b or ( II. 6b ) is simpler and suggests inter- 
esting interpretations, the natural framework for considering the variational problems i ll .41 
and jl.6\ is provided by time dependent families of Radon measures in R d . Following this 
point of view, one can replace p, by a continuous curve t 6 [0, 1] t—> p, (p, = p, in the 
absolutely continuous case) in the space M 1 J )C (R C ') of nonnegative Radon measures in 
endowed with the usual weak* topology induced by the duality with functions in 
The (Borel) vector field v, in \\.A\ induces a time dependent family of vector measures 
Vf := fi t \ t <C fit. In terms of the couple (p,V) the continuity equation d 1 -4b reads 

d t p, + V ■ V, = in the sense of distributions in &'(R d x (0, 1)), (1.12) 

and it is now a linear equation. Since v f = dv f /dp f is the density of V, w.r.t. fj.,, the action 
functional which has to be minimized in d 1 .4b can be written as 

p 

dp. (1.13) 



/ 


dv 


Jw 


dp 



Notice that in the case of absolutely continuous measures with respect to Jf d , i.e. p = p££ d 
and V = yvJ? d , the functional <t> p j can also be expressed as 



Pi i(jl,v) := / p .i(p,w)d^(x), p ,i(p,w):=p 

Jw 1 



(1.14) 



Denoting by C£(0, 1) the class of measure-valued distributional solutions p,V of the con- 
tinuity equation 01.12b , we end up with the equivalent characterization of the Kantorovich- 
Rubinstein- Wasserstein distance 

W/(po,pi) :=inf{^.i(M,v) : (p,v) e 66(0,1), M|, =0 = Mo, M|, =1 =Mi}- (1.15) 



Structural properties and convexity issues. The density function (j) = ty p \ : (0, +°°) x R d — > 
R d appearing in < | 1 - 1 4b exhibits some crucial features 

1. w i— » 0(-,w) is symmetric, positive (when w ^ 0), and p-homogeneous with respect 
to the vector variable w: this ensures that W p is symmetric and satisfies the triangular 
inequality. 

2. is jointly convex in (0, +°°) x M. d : this ensures that the functional 4> p .i (and therefore 
also S) defined in i ll. 13b is lower semicontinuous with respect to the weak* convergence 
of Radon measures. It is then possible to show that the infimum in d!.15b is attained, as 
soon as it is finite (i.e. when there exists at least one curve (p,v) 6 C£(0, 1) with finite 
energy S{fl,V) joining po to pi); in particular W^pojPi) = yields po = \l\. Moreover, 
the distance map (po,pi) i— > Wp(po,Pi) is lower semicontinuous with respect to the 
weak* convergence, a crucial property in many variational problems involving W p , as 

OH- 

3. (j) is jointly positively 1 -homogeneous: this a distinguished feature of the Wasserstein 
case, which shows that the functional <J> p . i depends only on p , v and not on the Lebesgue 
measure Jz? d , even if it can be represented as in 1 11. 14b . In other words, suppose that 
p = py and V = wy, where y is another reference (Radon, nonnegative) measure in M. d . 
Then 

<J> p .i(p,v) = / 4> p ,i(p,w)dy. (1.16) 

Jw 1 
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As we will show in this paper, the 1 -homogeneity assumption yields also two "quantita- 
tive" properties: if Ho is a probability measure, then any solution (jll,v) of the continuity 
equation ( 11.12b with finite energy S{ji,v) < +°° still preserves the mass p.t(M. d ) = 1 
for every time t > (and it is therefore equivalent to assume this condition in the def- 
inition of C£(0, 1), see e.g. (3.. Chap. 8]). Moreover, if the p-moment of /Xo m p (lJo) '■= 
J w i |jc| p dlio(x) is finite, then W p (p.o,fJ.i) < +°° if and only if m p (fj. l ) < +°°. 

Main definitions. Starting from the above remarks, it is then natural to consider the more 
general case when the density functional (j> : (0,+°°) x M. d — > [0,+°°) still satisfies 1. (p- 
homogeneity w.r.t. w) and 2. (convexity), but not 3. (1 -homogeneity). Due to this last choice, 
the associated integral functional <J> is no more independent of a reference measure y and it 
seems therefore too restrictive to consider only the case of the Lebesgue measure y = ££ d , 
In the present paper we will thus introduce a further nonnegative reference Radon mea- 
sure ye M 1 J )c (R rf ) and a general convex functional : (0,+°°) x W 1 — > [0,+°°) which is 
p-homogeneous w.r.t. its second (vector) variable and non degenerate (i.e. 0(p,w) > if 
w 0). Particularly interesting examples of density functionals <j>, corresponding to i ll. 6b , 
are given by 

IP 

(1.17) 

where h : (0, +°°) — > (0, +°°) is an increasing and concave function; the concavity of h 
is a necessary and sufficient condition for the convexity of <j> in J 1 . 17b (see [29 \ and Sj3]l. 
Choosing h(p) := p™, a 6 (0, 1), one obtains 

p 

= p e -P\yv\ p , d:=(l-a)p + ae(l,p), (1.18) 

which is jointly ^-homogeneous in (p,w). 

In the case, e.g., when a < 1 in dl.l8b or more generally lim p | 00 /!(p)/p = 0, the reces- 
sion function of (j) satisfies 

f°(p,w)= lim X- 1 <j>(Xp,Xvi)=+oc ifp.w^O, (1.19) 

AT+°° 

so that the associated integral functional reads as 

4>0,v|y) := / 0(p,w)dy i u=P7+M ± 5 v = wy<y, (1.20) 

extended to +°° when V is not absolutely continuous with respect to y or supp(ii) <f_ supp(y). 
Notice that only the density p of the y-absolutely continuous part of /i enters in the func- 
tional, but the functional could be finite even if [I has a singular part J u ± . This choice is 
crucial in order to obtain a lower semicontinuous functional w.r.t. weak* convergence of 
measures. The associated (0,y)-Wasserstein distance is therefore 

W; r ( Mo ,Mi):=inf{^,r(M,v):( i u,v)ee£(0,l), M|, =0 = Mo, M| f=1 =Mi}, (1-21) 
where the energy S^ y of a curve (/i, v) e C£(0, 1) is 

4. 7 (M,v):= [ 0(^,V t \y)6t. (1.22) 



4>(p,w) :=%) 



h(p) 



w 
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The most important case associated to the functional ( 11.18b deserves the distinguished nota- 
tion 

The limiting case a = 8 = 1 corresponds to the ZZ -Wasserstein distance, the Sobolev Wy X ' p 
corresponds to a = 0, 9 = p. The choice of y allows for a great flexibility: besides the 
Lebesgue measure in W 1 , we quote 

- 7 := \ n , ■Q being an open subset of R d . The measures are then supported in Q, and, 
with the choice (1.171 and v = y//h{p), d 1 - 12b is a weak formulation of the continuity 
equation (n^Q being the exterior unit normal to dQ) 

d,pt+V-(h(p t )v t ) =0 in 12 x (0,1), Vrn 3fi =0 on (9X2. (1.24) 

This choice is useful for studying equations dl.9l l (see II II ). 1 11. 10b , 1 11. 1 lb in bounded 
domains with Neumann boundary conditions. 

- 7 := e- y ^ d for some C 1 potential V : R d -> R. With the choice (THl and v = w/h(p) 
J 1 . 1 2b is a weak formulation of the equation 

dtPt+V-(h(ftt)vt)-h(pt)W-v t = inR rf x (0,1). (1.25) 

When h(p) = p a , p = 2, the gradient flow of := ^-g^i-g) f 2 " is ^ 

Kolmogorov-Fokker-Planck equations 1 15] 

d,H -An - V ■ (mVV) = 0, a t p - Zip + VV ■ Vp = 0, 

which in the Wasserstein framework is generated by the logarithmic entropy ( I19l l3ll5l). 

- 7 := -3V k \,„, M being a smooth ^-dimensional manifold embedded in R d with the Rie- 

' |M ° 

mannian metric induced by the Euclidean distance; J^ k denotes the A^-dimensional Haus- 
dorff measure, i ll. 12b is a weak formulation of 

d t p, +div M (/*(p)v,) =0 onMx(0,l). (1.26) 

Thanks to Nash embedding theorems 1221231 . the study of the continuity equation and of 
the weighted Wasserstein distances on arbitrary Riemannian manifolds can be reduced 
to this case, which could be therefore applied to study equations i ll. 9b . dl.lOb . dl.llb on 
Riemannian manifolds. 

Main results. Let us now summarize some of the main properties of W Pt a;y(- , *) we will 
prove in the last section of the present paper. In order to deal with distances (instead of 
pseudo-distances, possibly assuming the value +°°), for a nonnegative Radon measure a 
we will denote by M p . a; y[o] the set of all measures p with Wp, a; y(p,CT) < +°° endowed 
with the Wyj.c^y-distance. 

1- Mj^yja] is a complete metric space (Theorem l5.7b . 

2- Wp <a -y induces a stronger convergence than the usual weak* one (Theorem [53}. 

3. Bounded sets in M p .a;y[<7] are weakly* relatively compact (Theorem l5.5b . 

4. The map (po,pi) <— > Wp jCt; y(po,pi) is weakly* lower semicontinuous (Theorem 15.6b . 
convex (Theorem 1 5. lib , and subadditive (Theorem 15. 12b . It enjoys some useful mono- 
tonicity properties with respect to y (Proposition 15.14b and to convolution (Theorem 

E33. 
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5. The infimum in dl.l5l > is attained, M^.ce ;r [a] is a geodesic space (Theorem 15.41 ), and 
constant speed geodesies connecting two measures [1q,[1i are unique (Theorem 15.1 lb . 

6. If 

/ W-^ e - 1 )d 7 W<+oo e = (i-a) P +a, -£- = -i-, (1.27) 

iW>i 0-1 1-a 

and ae?^), then M p , a;y[c] C 3 > (]R C ') (Theorem l5.8b . If moreover / satisfies stronger 
summability assumptions, then the distances W p , a ^ provide a control of various mo- 
ments of the measures (Theorem 15. 9b . Comparison results with W p and W ' F are also 
discussed in £15.41 

7. Absolutely continuous curves w.r.t. W p . a; y can be characterized in completely analogous 
ways as in the Wasserstein case ( i]5.3b . 

8. In the case y = 5£ d the functional 

%(,\y) := (2 _ a ; (1 _ a) / Rrf P 2 -"dx M = 9* « (1.28) 

is geodesically convex w.r.t. the distance W 2 a .^,i and the heat equation in M. d is its 
gradient flow, as formally suggested by dl.9b ( 35.51 we prove this property in the case 
a > 1 - 2/d, when ^(W 1 ) is complete w.r.t. W 2 a .&i-) 

Plan of the paper. Section 2 recalls some basic notation and preliminary facts about weak* 
convergence and integral functionals of Radon measures; [23] recalls a simple duality result 
in convex analysis, which plays a crucial role in the analysis of the integrand 0(p,w). 

The third section is devoted to the class of admissible action integral functionals <P like 
i ll. 20b and their density (j). Starting from a few basic structural assumptions on (j) we deduce 
its main properties and we present some important examples in Section[3]2] The correspond- 
ing properties of 4> (in particular, lower semicontinuity and relaxation with respect to weak* 
convergence, monotonicity, etc) are considered in Section [3~3l 

Section|4]is devoted to the study of measure- valued solutions of the continuity equation 
i ll. 12b . It starts with some preliminary basic results, which extend the theory presented in 
|3i] to the case of general Radon measures: this extension is motivated by the fact that the 
class of probability measures (and therefore with finite mass) is too restrictive to study the 
distances Wp.a-y, in particular when y(R rf ) = +°° as in the case of the Lebesgue measure. 
We shall see (Remark |5.27b that 9(1^) with the distance W p a .^>,i is not complete if d > 
p/{Q — I) = q/{\ — a). We consider in Section l4~2l the class of solutions of i ll. 12b with finite 
energy l ll.22b . deriving all basic estimates to control their mass and momentum. 

As we briefly showed, Section [5] contains all main results of the paper concerning the 
modified Wasserstein distances. 



2 Notation and preliminaries 

Here is a list of the main notation used throughout the paper: 

Br The open ball (in some M. ) of radius R centered at 

2(R ; ') (resp. B C (R A )) Borel subsets of R h (resp. with compact closure) 

5 > (R / ') Borel probability measures in R /! 

M^IR' 5 ) (resp. M^M' 5 )) Finite (resp. Radon), nonnegative Borel measures on K 

y(R /! ) Borel probability measures in K 
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M(R /, ;R m ) R m -valued Borel measures with finite variation 

Mioc (K A ; R m ) R m -valued Radon measures 

1 1 jU 1 1 Total variation of ju E M ioc (R A ; W" ) , see d2T2l 

C®(M. h ) Continuous and bounded real functions 

m p ( J u) p-moment f Rrf \x\ p dp of jX e M + (M /l ) 

y/°° Recession function of \ff, see l |2.4| l 

f(ju|7), <&(/i,v|y) Integral functionals on measures, see l2.2l and l3.3l 

(/*,£) the integrals VCdM,VC-djU 

C£(0,r), ££,)). y(0,T), Classes of measure-valued solutions of the continuity 
£8.(0, T;hq jui) equation, see Def.|43]and Sec.l4~2l 



2.1 Measures and weak convergence 

We recall some basic notation and properties of weak convergence of (vector) radon mea- 
sures (see e.g. (2J). A Radon vector measure in M loc {R h ;M. m ) is a Revalued map ju : 
S c (R /! ) — > R m defined on the Borel sets of R' 1 with compact closure. We identify ju 6 
Mi oc (R /! ;R m ) with a vector (ju 1 ,ju 2 , ■ ■ • ,/i m ) of m measures in Mi oc (R /l ): its integral with a 
continuous vector valued function with compact support £ 6 C^(R ;R m ) is given by 

<M,C}:= / C ■ dM = £ ! C\x)d l i i {x). (2.1) 

It is well known that Mi oc (R /! ;R m ) can be identified with the dual of C°(R ft ;R m ) by the 
above duality pairing and it is therefore endowed with the corresponding of weak* topology. 
If || • || is a norm in R'' with dual || • ||* (in particular the euclidean norm • |) for every open 
subset A C R A we have 

\\H\\{A) = mv{j^.dn: supp(£)cA, ||C(x)||*<l VieR"}. (2.2) 

\\fl\\ is in fact a Radon positive measure in M 1 J )c (R /l ) and ju admits the polar decomposition 
/I = w\\fl\\ where the Borel vector field w belongs to Li 0C (|[/i||;R m ). We thus have 

<fi,C>= / C-d M = / C-wd||/i||. (2.3) 

JR* jR h 

If (jU /t ) / tgN is a sequence in Mi oc (R /, ;R m ) with sup n ||ju||(5«) < +°° for every open ball 
Br, then it is possible to extract a subsequence fit weakly* convergent to ju 6 M(R /l ;R m ), 
whose total variation \\n kn \\ weakly* converges to X 6 M + (R A ) with ||ju|| < A. 



2.2 Convex functionals defined on Radon measures 

Let \j/ : R m — > [0, +°°] be a convex and lower semicontinuous function with \jf{0) = 0, whose 
proper domain D(\j/) := {x 6 R m : < +°°} has non empty interior. Its recession function 
(see e.g. (2)) : R m -> [0, +°°] is defined as 

rW:= lim V^Zl =sup ZM. (2 .4) 
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\ff°° is still convex, lower semicontinuous, and positively 1 -homogeneous, so that its proper 
domain Z)(ty/°°) is a convex cone always containing 0. We say that 

\ff has a superlinear growth if y°°(y) = °° for every y ^ 0: D(\ff°°) — {0}, 
\j/ has a sublinear growth if y°°(y) = for every y 6 K m . 

Letnow ye M 1 + c (]R /l ) and/i £ Mi oc (]R /l ;]R m ) withsupp(ju) C supp(y); the Lebesgue decom- 
position of fi w.r.t. y reads ju = #y+jU x , where # = dju/dy. We can introduce a nonnegative 
Radon measure CT € M 1 + c (M /! ) such that /i x = ■& ± a <S cr, e.g. a = |/X X | and we set 

W a (jl\y) := [ w(&(x))dy(x), f~(/l|y) ■= [, Vr{* x (y))d0<y), (2.6) 

JR* JR h 

and finally 

•P(Mlr) :=f a (/i|y)+ , P°°(M|y); f°°0|y)=+°°if supp(ju) £ supp(y). (2.7) 

Since is 1 -homogeneous, the definition of f°° depends on y only through its support and 
it is independent of the particular choice of a in d2.6b . When i|/ has a superlinear growth then 
the functional f is finite iff ju <C y and f fl (/i|y) is finite; in this case f(^|y) = f fl (ju|y). 

Theorem 2.1 (L.s.c. and relaxation of integral functionals of measures 1 1] 2|) Let us con- 
sider two sequences y„ e Mj+ C (R A ) , /i„ e Mi oc ; R m ) weakly* converging to y e M jj^ (M. h ) 
and /i 6 Mi oc (R A ;R'") respectively. We have 

Hminf'P(/i n |y n )>'P(M|y). (2.8) 

«T+°° 

Lef conversely fl,ybe such that f (ju|y) < +°°. 77zen rtere exists a sequence jU n = #„y <C y 
weaWy* converging to fl such that 

lim f fl ( M „ | y) = lim / (*) ) dy(x) = | y) . (2.9) 

Theorem 2.2 (Montonicity w.r.t. y)Ify\ < fi then 

f(/i|y2)<f(M|yi)- (2.10) 

Proof Thanks to Theorem |2.H it is sufficient to prove the above inequality for ju <g; y 1 . Since 
y = Qyj, with density 6 < 1 yi-a.e., we have /X = D'f with ■& 1 = d and therefore 

/ iK*i)dyi=/ v(e- l &2)edy 2 > [ wi^dyz, C2.ii) 

where we used the property ©^(f^'jc) > for < 1, being )|/(0) =0. □ 

Theorem 2.3 (Monotonicity with respect to convolution) Ifk e C"(K C ') is a convolution 
kernel satisfying k(x) > 0, f R( / k(x) dx = 1 , f/;en 

f(/i*/fc|y*£) <f(ju|y). (2.12) 

The proof follows the same argument of [ 3 ; , Lemma 8.1.10], by observing that the map 
(x,y) i — > x\jf(y/x) is convex and positively 1-homogeneous in (0, +°°) x W 1 . 
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2.3 A duality result in convex analysis 

Let X, Y be Banach spaces and let A be an open convex subset of X. We consider a convex 
(and a fortiori continuous) function <j> : A x Y — » R and its partial Legendre transform 

^(x,y*):=sup<y*,y)-^(*,y)e(-«»,+oo], VjteA./er*. (2.13) 

The following duality result is well known in the framework of minimax problems [29 \. 
Theorem 2.4 is a l.s.c. function and there exists a convex set Y* C F* swcn fnaf 

ft*,/) <-h» «■ /ey;, (2.U) 

so f/iaf ^(•,3'*) = +°°for every y* 6 T* \ T * and admits the dual representation formula 
4>(x,y) = sup (y,/)-£(xy) VxeA,yeT. (2.15) 

For every y* 6 F * we nave 

the map xt-^<p(x,y*) is concave (and continuous) in A. (2. 16) 

Conversely, a function (j) : A x 7 — > R /'.? convex j/if admits the dual representation d2. 15b for 
a function (j) satisfying d2.16l l. 

Proof Let us first show that d2.16l l holds. For a fixed y* e Y*, xq,x\ > 0, 9 6 [0, 1], and 
arbitrary 6 7, we get 

0((i-0)* o + 0*i,y*)> </,(i-t%o + #yi> + -i?)y + i?yi) 
>(l-^)((/,yo)-<?(x ,yo))+^((/,yi)-<?'(^,yi)). 

Taking the supremum with respect to yo,yi we eventually get 

<K(1 - &)x + Vx h y*) > (1 - #)0(*o,/) + #<M*i,y*) (2.17) 

and we conclude that <j>(-,y*) is concave. In particular, if it takes the value +°° at some point 
it should be identically +°° so that d2. 14b holds. 

The converse implication is even easier, since d2. 15b exhibits (j) as a supremum of con- 
tinuous and convex functions (jointly in x e A,y 6 Y). □ 

3 Action functionals 

The aim of this section is to study some property of integral functionals of the type 

4> a (M,v|7) := / <Hp,w)dy, fi = pye M+ c (R d ), v = wye M loc (R d ;R d ) (3.1) 
and their relaxation, when <j> satisfies suitable convexity and homogeneity properties. 
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3.1 Action density functions 

Let us therefore consider a nonnegative density function (j) : (0, +°°) x R — > [0, +°°) and an 
exponent p 6 ( 1 , +°°) satisfying the following assumptions 

(j) is convex and (a fortiori) continuous, (3.2a) 

w I— > w) is homogeneous of degree /?, i.e. 

, (3.2b) 
0(p,Aw) = |A|^(p,w) Vp >0,A eR,weR rf , 

3po>0: 0(por) is non degenerate, i.e. 0(p o ,w)>O Vw6R £ '\{0}. (3.2c) 

Let q = p/{p— 1) 6 (1,+°°) be the usual conjugate exponent of p. We denote by : 
(0, +°°) xl^-t (— °°, +°°] the partial Legendre transform 

~0(p,z) := sup z-w-~0(p,w) Vp>0,zeR J . (3.2d) 
# weR'' P 

We collect some useful properties of such functions in the following result. 
Theorem 3.1 Let (j) : (0, +°°) x R d -> W 1 satisfy &3k.b.c). Then 

1. For every p > the function w i— > (j) (p , w) l / p is a norm ofW 1 whose dual norm is given 
byz^ ^(PjZ) 1 /*, i.e. 

^(p i z) 1 ^ = sup W ' Z ^(p,w)^ = SU p J"' (3.3) 

w^o ^(p^) 1 //' Z7 , 0(p,z)V« 

/« particular <p(-,z) ls q-homogeneous with respect to z. 

2. 77;e marginal conjugate function fafe.? i?s values in [0, +°°) and for every z 6 R rf 

rfte map p i— > 0(p,z) is concave and non decreasing in (0, +°°). (3.4) 
/« particular, for every w £ R rf 

f/;e map p i— > (p , w) is convex and non increasing in (0, +°°) . (3.5) 



0(p, z) < (a + fcp) |z| ? , <^)(p, z) > (a + fcp) 1 p |w| p Vp>0,z,weR rf . (3.6) 

4. For every closed interval [po , pi ] C (0, +°°) //iere exists a constant C = Cn,,^ > such 
that for every p 6 [po,Pl] 

C^Vl* < 0(p,w) <C|wj p , C~ l \z\<i <${p,z) <C\z\ q Vw,z£R d . (3.7) 

Equivalently, a function (j) satisfies d3.2h .b.c) if and only if it admits the dual representation 
formula 

-(j>(p,Vf) = sup wz- ~0(p,z) Vp>0,weK'', (3.8) 



where <j) : (0, +°°) x R d — > (0, +°°) is a nonnegative function which is convex and q-homogeneous 
w.r.t. z and concave with respect to p. 
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Proof Let us first assume that <j) satisfies d3.2b ,b,c). The function w i— » 0(p,w) 1 / p is 1- 
homogeneous and its sublevels are convex, i.e. it is the gauge function of a (symmetric) 
convex set and therefore it is a (semi)-norm. The concavity of follows from Theorem 
12.41 taking w = in d3.2dl i. we easily get that is nonnegative; d3.2cb yields, for a suitable 
constant cq > 0, 

<KPo,w) >c |w| p Vwel 11 , so that 0(p o ,z) < c |z| 9 < VzeR rf . (3.9) 

Still applying Theorem [2.4l we obtain that p i— > 0(p,z) is finite, strictly positive and nonde- 
creasing in the interval (0, +°°). Since 0(p,O) = we easily get 

^(p,z)<^(po,z)< Co |z|« VzeR'', pe(0,p ); (3.10) 

§(P,z)<- ^(Po,z)< — p|z|« Vz€l J , pe(p ,+°°). (3.11) 
Po Po 

Combining the last two bounds we get d3.6l l. d3.7b follows by homogeneity and by the fact 
that the continuous map (j) has a maximum and a strictly positive minimum on the compact 
set [p ,pi] x{wel lf : |w| = 1}. 

The final assertion concerning d3-8b still follows by Theorem l2.4l □ 



3.2 Examples 

Example 3.2 Our main example is provided by the function 

2 

02,a(p,w) = -L-L fa a (p,z):=p a \z\ 2 , 0<a<l. (3.12) 
p" 

Observe that fe.a is positively 0-homogeneous, d := 2 — a, i.e. 

^2.a(Ap,Aw) = A e 0(p,w) VA,p >0, weMf 1 . (3.13) 
It can be considered as a family of interpolating densities between the case a = 0, when 

<fc,o(p,w) := |w| 2 , (3.14) 
and a = 1, corresponding to the 1-homogeneous functional 

2 

<te,i(P,w):=LL. (3.15) 

Example 3.3 More generally, we introduce a concave function h : (0, +°°) — > (0, +°°), which 
is a fortiori continuous and nondecreasing, and we consider the density function 

2 

0(p,w) := i^y, 0(p,*):=A(p)|w| 2 . (3.16) 

If & is of class C 2 , we can express the concavity condition in terms of the function g(p) := 
l/h(p) as 

A is concave & g"(p)g(p)>2(g'(p)) 2 Vp > 0, (3.17) 

which is related to a condition introduced in |6 Section 2.2, (2.12c)] to study entropy func- 
tionals. 
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Example 3.4 We consider matrix-valued functions H, G : (0, +°°) — > M. dxd such that 

H(p),G(p) are symmetric and positive definite, H(p) = G _1 (p) Vp>0. (3.18) 
They induce the action density (j) : (0, +°°) x K ' — > [0, +°°) defined as 

<Hp,w) := (G(p)w,w> = <H- 1 (p)w,w) . (3.19) 
Taking into account Theorem [XT] <j> satisfies conditions 0.21 ) if and only if the maps 

p i — > (H(p)w,w) are concave in (0,+-) VweK d . (3.20) 

Equivalently, 

H((l-#)p + tfpi) > (l-iJ)H(po) + iJH(p 1 ) as quadratic forms. (3.21) 

When G is of class C 2 this is also equivalent to ask that 

G"(p)>2G'(p)H(p)G'(p) Vp>0, (3.22) 

in the sense of the associated quadratic forms. In fact, differentiating H = G~' with respect 
to p we get 

H' = -HG'H, H" = -HG"H + 2HG'HG'H, 

so that 

d 2 

-p-^ (H(p)w,w) = — (G"w,w) +2(G'HG'w, w) where w := Hw; 
we eventually recall that H(p) is invertible for every p > 0. 

Example 3.5 Let || ■ || be any norm in W 1 with dual norm || ■ ||*, and let h : (0, +°°) — > (0, +°°) 



be a concave (continuous, nondecreasing) function as in Example l3.3l We can thus consider 

Hp,z):=h(p)\\z\\l (3.23) 



4>(p,w) :=h(p) 



w 



h(p) 

See [20 21 1 for a in-depth study of this class of functions. 

Example 3.6 ((a-O)-homogeneoits junctionals) In the particular case h(p) := p a the func- 
tional (j) of the previous example is jointly positively 0-homogeneous, with 8 := a + (1 — 
Oc)p. This is in fact the most general example of ^-homogeneous functional, since if (j) is 
^-positively homogeneous, 1 < 6 < p, then 

<Kp,w)=p 9 <Kl,w/p)=p e -^(l,w)=p a ||w/p a |r, a= EZ ^-, (3-24) 

p-l 

where ||w|| := 0(l,w) 1//p is a norm in M. d by Theorem 13. II The dual marginal density ^ in 
this case takes the form 

$(p,z)=p a \\z\\l Vp>0, zeW', (3.25) 
and it is q + a-homogeneous. Notice that a and 9 are related by 

- + - = 1. (3.26) 

p q 

In the particular case when || ■ || = || ■ ||„ = | ■ | is the Euclidean norm, we set as in d3.16l > 



w 



fp.a 



(p,w):=p« — 



V(p,z):=p a |z| 9 , 0<a<l. (3.27) 
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3.3 The action functional on measures 

Lower semicontinuity envelope and recession function. Thanks to the monotonicity prop- 
erty d 3 - 5 b . we can extend (j) also for p = by setting for every w € M. d 

<MO,w)=sup<Kp,w)=lim<Hp,w); in particular J 0) " °' (3.28) 
p>0 PlO ^<p(0,w)>0 rfw^O. 

When p < we simply set 0(p,w) = +°°, observing that this extension is lower semicon- 
tinuous in R x R rf . It is not difficult to check that 0(0, ■) satisfies an analogous formula 

0(O,z) = sup z-w-0(O,w) = inf 0(p,z) = lim0(p,z) VzeR rf . (3.29) 

weR'' P >0 ^° 



Observe that, as in the (a-0)-homogeneous case of Example 1 3. 6 1 with a > 0, 

c>;0.z; = - l() . w ,-J+~ * W ^° (3.30) 

10 ifw = 0. 

As in | |2.4| |. we also introduce the recession functional 

0°°(p.w) = sup ^0(Ap,Aw) = lim i(^(Ap,Xw)= lim A^'^Ap.w). (3.31) 

0°° is still convex, p-homogeneous w.r.t. w, and l.s.c. with values in [0, +<»]; moreover, it is 
1 -homogeneous so that it can be expressed as 

(P.w) = < p/ (3.32) 
[ +°° if p = and w ^ 0, 

where (p°° : R d — > [0, +°°] is a convex and p-homogeneous function which is non degenerate, 
i.e. <p°°(w) > if w ^ 0. <p°° admits a dual representation, based on 

9~(z) := inf = lim \§{Xp,z). (3.33) 

9 00 is finite, convex, nonnegative, and g-homogeneous, so that <p°°(z) l ' q is a seminorm, 
which does not vanish at z £ R d if and only if p i— > (p , z) has a linear growth when p | +°°. 
It is easy to check that 

<p°°(yv) 1/p = sup jwz : <p°°(z) < l|. (3.34) 

In the case (j) has a sublinear growth w.r.t. p, as for (a-0)-homogeneous functionals with 
a < 1 (see Example l3.6b , we have in particular 

f+oo if w =^0 
when ?»=<), ?"(w) = ^ ^ ' (3.35) 

if w = 0. 
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The action functional. Let e M 1 J )c (IR rf ) be nonnegative Radon measures and let V 6 
Mi 0C (R c ';R'') be a vector Radon measure onR . We assume that supp(/i),supp(v) C supp(y), 
and we write their Lebesgue decomposition with respect to the reference measure 7 

M:=py + i u ± , v:=wy+v x . (3.36) 

We can always introduce a nonnegative Radon measure a £ M + (I2) such that J u ± = p X (7 <C 
cr,V x = w x o~ <C O, e.g. <7 := jU x + V x |. We can thus define the action functional 

4>( J u,v|7) = $ fl (M,v|7) + 4>~( J u,v|7):= / *(p,w)dy+ / 0~(p ± ,w ± )da. (3.37) 

Observe that, being 0°° 1 -homogeneous, this definition is independent of o~. We will also 
use a localized version of 4>: if B 6 25 (R d ) we set 

*(M,V|y,B) := / 0(p,w)dy+ / oc (p ± ,w ± )da. (3.38) 



Lemma 3.7 Lef ,u = py+ i u ± , V = wy+ V x be such that <5(ji,v|y) is finite. Then V x = 
w x ,u x <C jU x a«rf 

*-( M ,v|y)= / (< p~(w x )d i u x , *(M,v|y)= / *(p,w)dy+ / ; < p"(w x )d i u x . (3.39) 

Moreover, if <j> has a sublinear growth with respect to p (e.g. in the (a-d)-homogeneous 
case of Example \3.6\ with a < 1) then (f>°°(-) = and 

#(M,V) < +°° V = w 7< 7, #(/i,V) = * a O,v) = / 0(p,w)dy, (3.40) 

independently on the singular part fl^. 

Proof Let ct 6 M 1 "J K .(R </ ) any measure such that jU x <C ct, |v x <S CT so that 4 , ° ( J u,v|y) can 
be represented as 

*-(M,v|r) = / rj ^(p x ,w x )da, p x = ^, w x = ^. 

When <J , °° (/i , V | y) < +°°, ( 13.32b yields w x (x) = for CT-a.e. x such that p x (x) = 0. It follows 
that 

<P(H,v)<+°° => V x </rS (3-41) 

so that one can always choose a = /x x , p x = 1, and decompose V x as w x /x x obtaining 
d3.39b - J3.40l > is then an immediate consequence of d3.35b . □ 

Remark 3.8 When 0(0, z) = (e.g. in the (a-0)-homogeneous case of Example 13.61 with 
a > 0) the density w of V w.r.t. y vanishes if p vanishes, i.e. 

<P(H,v\y) <+°° => w(jc) =0ifp(jc) =0, for y-a.e. xeM. d . (3.42) 

In particular V a is absolutely continuous also with respect to pL. 

Applying Theorem [2T| we immediately get 
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Lemma 3.9 (Lower semicontinuity and approximation of the action functional) The 

action functional is lower semicontinuous with respect to weak* convergence of measures, 
i.e. if 

p„^*p, 7,^*7 weakly* in M,+ c (R d ), V„^*V in Mi oc (K d ;H d ) as n | +°°, 
then 

limM<&(p„,v„|y„) > #(/i,v|y). 

Equiintegrability estimate. We collect in this section some basic estimates on <j> which will 
turn to be useful in the sequel. Let us first introduce the notation 

\\z\\*--$(l,z) l /<i, \\yr\\:=4>(l,vr)VP, ij-'M < ||z||, < tj|z|, (3.43) 

L} := \(a,b) : sup 0(p,z) < a + bpX, h(p) := inf la + bp : (a, b) e L}\, (3.44) 

IMI»=i 

H(s,p) :=sh{p/s) = inf \as + bp : [a,b)&rA. (3.45) 

Observe that h is a concave increasing function defined in [0, +°°), satisfying, in the homo- 
geneous case h(p) = h(p) = p a . It provides the bounds 

£(p,z)<h(p)||z||f, ||w||<h(p) 1 /^(p,w) 1 ^, 

9 00 (z)<h 00 ||z||«, || w|| < (h-) 1/? 9"(w) 1 /", if h~ := lim A -1 h(A) > 0. (3A6) 

Observe that when h°° = then <p°° = and <p°°(w) is given by ( 13.35b . 

Proposition 3.10 (Integrability estimates) Let £ be a nonnegative Borel function such that 

:= / C 9 dp and y(C 9 ) := / C 9 dy arefinite, 

J«. d Jm. d 

andletZ := {x e R d : £(jc) >0}. //<£(p,v|y) < +°° we have 

, CWd||v||W < ^ 1 /P( M)V |y ) Z)H 1 /?( 7 (^) )jU (^)). (3.47) 

/« particular, for every Borel set A 6 33 (R rf ) we /zave 

||v||(A) < 4> 1 ^(p,v|y,A)/f 1/9 (y(A),p(A)). (3.48) 

Prao/ It is sufficient to prove d3.47b - Observe that if (a,b) 6 rj, then a > 0, and h°° < b so 
that by d3.46l > we have 



CWd|jv||W< / C|M|dy+ / CIKHdp 

rf JZ JZ 

< (^(p,w)dy) 1/ "(^C 9 h(P)dy) lA ' + ( j (<p 0O (w ± )d M ± ) 1/ ''(h^C 9 dM 1 " 1 " 



< (<f (p,v|y,Z)) P (aJ w ? A y +b J^) 
Taking the infimum of the last term over all the couples (a,&) £ i~0 we obtain ( 13.48b . □ 
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4 Measure valued solutions of the continuity equation in Mr 

In this section we collect some results on the continuity equation 

d,fi,+V-v,=0 inR d x(0,r), (4.1) 

which we will need in the sequel. Here [1,,V, are Borel families of measures (see e.g. |3]) in 
Mj+ c (R rf ) and Mi oc (R d ;M'') respectively, defined for t in the open interval (0, T), such that 



M,(5«)df<+~ V R := \v,\(B R )dt <+™ VR > 0, (4.2) 

Jo 

and we suppose that fl4. lb holds in the sense of distributions, i.e. 



d,£(x,t)dn,{x)dt+ / V,C(x,0-dv,(x)d» = (4.3) 

for every C e (R* 7 x (0, J)). Thanks to the disintegration theorem 031 4, 111-70], we can 
identify (V() J6 m,r) w i m me measure V = / r v,dr e Mi oc (R rf x (0,T);M d ) defined by the 
formula 

(V,Q = £ (^(x^-dv,^ dt VCGC?(ffi d x(0,r);R d ). (4.4) 



4.1 Preliminaries 

Let us first adapt the results of [3 Chap. 8] (concerning a family of probability measures 
fi t ) to the more general case of Radon measures. First of all we recall some (technical) 
preliminaries. 

Lemma 4.1 (Continuous representative) Let pt, , V, be Borel families of measures satisfy- 
ing d4.2| | and 04.3I I. Then there exists a unique weakly* continuous curve t g [0, T] i— > jl, 6 
Mj^ (R d ) sucn r«a? il t = p,, for J? l -a.e.te (0, T);ifC,e C x c {M d x [0, T] ) ant/ f i < t 2 € [0, T], 
we have 

[ Ct 2 d £, 2 - / C, dfti = f' 2 [ d,t; dm (x) df + f / ■ dv, ( x ) dt, (4.5) 

and f/ie majs of p., can be uniformly bounded by 

sup fit (Br) < p s {B 2R )+2R- l V 2R \/s e [0,T]. (4.6) 

te[o,T] 

Moreover, if fi. s (M. d ) < +°°for some s 6 [0, 7"] and \m\ R ^ +co R^ l V R = 0, then the total mass 
p,t(M d ) is (finite and) constant. 

Proof Let us take f (x,f) = ?](*)£(*)> H € C~(0,r) and £ e C^°(R d ) with supp£ C 5«; we 
have 

" f C W dMr (*)) * = jf H(0 (/ Rrf VCW • dv f (x)) dt, 

so that the map 1 1— > /i r (0 = Jr^ Cd/x, belongs to W 1,1 (0, 7*) with distributional derivative 
MO=/ V;(x)-dV,(x) fijrJSr 1 -a.e.te(0,r), (4.7) 
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satisfying 

|M,(C)|<V R (f)sup|V£|, V R (t):=\v,\(B R ), [ V R (t)dt = V R < +<». (4.8) 

If is the set of its Lebesgue points, we know that «Jzf 1 ((0, T)\Lr) =0. Let us now take 
an increasing sequence R„ :=2" f +°° and countable sets Z„ C C™(B Rn ) which are dense in 
C^(B R J := {£ € C 1 (R </ ) : supp(£) C 5^"}, the closure of Cl(B Rn ) with respect the usual 
C 1 norm ||£|| c i = sup ffi( /(|£|, |V£|). We also set Lz '■= n„ e p^ 6Zn L£. The restriction of the 
curve ji to Lz provides a uniformly continuous family of functionals on each space Cq (B Rn ), 
since l !4.8l l shows 

\li t (Q-Li s (0\<K\yfv Rn (X)M Vj.teiz VCeZ„. 

Therefore, for every n 6 N it can be extended in a unique way to a continuous curve 
{£r"}/e[o.r] m [Co(5r (i )]' which is uniformly bounded and satisfies the compatibility con- 
dition 

ft"(C)=A"(f) ifm^nandCGQ 1 ^)- (4-9) 
If C € C\ [W 1 ) we can thus define 

:= jU,"(0 for every n e N such that supp(^) C 5fi„. (4.10) 

If we show that {jU, (S« n ) } (6 l z is uniformly bounded for every n 6 N, the extension provides 
a continuous curve in Mjt (K d ). To this aim, let us consider nonnegative, smooth functions 

: R rf -> [0, 1], such that := £b(*/2*), (4.1 la) 

^(jc) = 1 if |x| < 2 fe , &(x) = 0if |jc| >2* +1 , |V5tWI<A2 - *, (4.11b) 

for some constant A > 1 . It is not restrictive to suppose that £j £ Zj. + j . Applying the previous 
formula d4.7b . for f , s £ Lz we have 

IMSt) " ft(fit)| < «* := 2 ] - k £ |V,| (B2R t dr < A2~ k V 2Rk . (4.12) 
It follows that 

H,(B Rk ) < < Hs(Ck) +A2- k V 2Rk < Vs{B 2Rk ) +A2- k V 2Rk Vr £ L z . (4.13) 
Integrating with respect to j we end up with the uniform bound 

H,(B Rk ) < A2- k V Rk+1 +j\ s {B 2Rk )ds < +oo \/ t e L z . 

Observe that the extension p., satisfies d4.13l > (and therefore, in a completely analogous way, 
d4~6l >> and J4TT21 for every s , t £ [0, T] . 

Now we show d4.5l l. Let us choose C, £ Cj. (R d x [0, r] ) and 7j e £ C" (fi , t 2 ) such that 

0<t/ 8 (0<1, HmT]e(t)=X(, l<h )(t) Vre[0,J], lim^ = ^-^ 
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in the duality with continuous functions in [0, T] . We get 

0= / / d f (T} £ £)d J u f (x)df+ / / V^Q-dVtdt 

Jo JW' Jo Jr j 

= [ T) £ (f) / <9,£d/x,dr+ / rj e (t) f V x C-dv,dt+ f n' e {t) f CdMrdf. 
Jo Jw 1 Jo Jw 1 Jo Jk 1 

Passing to the limit as e vanishes and invoking the continuity of p.,, we get d4.5l l. 

Finally, if lim R f +00 W _1 Vr = we can pass to the limit as f +°° in the inequality 
(14. 12b . which also holds for every t,s E [0, T] if we replace fi by fx, by choosing s so that 

m := £ s (R f ') = lim < +oo. 

It follows that fi,(M. d ) = lim iT+OOJ Qi(^) = m for every f e [0, T), □ 



Thanks to Lemma 14.11 we can introduce the following class of solutions of the continuity 
equation. 

Definition 4.2 (Solutions of the continuity equation) We denote by QE(0,T) the set of 
time dependent measures (Ht)te[0,T]j i v t)te{0,T) sucn mat 

1. t i — ► jtXf is weakly* continuous in M 1 J )C (K £/ ) (in particular, sup fg j r]M/(^«) < +°° f° r 
every R > 0), 

2. (v ( ) te ( 0l r) is a Borel family with J \v, \(B R )dt < +<=° VR > 0; 

3. (jU,V) is a distributional solution of ( 14.1b . 

C£ (0, T,o — ► r\) denotes the subset of (ii,v) e S£(0, J) such that iio = ct, jUi = r). 
Solutions of the continuity equation can be rescaled in time: 

Lemma 4.3 (Time rescaling) Let t : s £ [0, r'] — > t(s) e [0, T) be a strictly increasing ab- 
solutely continuous map with absolutely continuous inverse s := t . Then [jl, v) is a distri- 
butional solution of (14. 1 b if and only if 

p := /I o t, V := t (v o t) , is a distributional solution of d4- 1 b on (0, 7" ) . 

We refer to [3 , Lemma 8.1.3] for the proof. 

The proof of the next lemma follows directly from l(4 



Lemma 4.4 (Glueing solutions) Let (ji', V ; ) e C£ (0, 7}), i = 1 , 2, wif/t ^ = jUq. Then the 
new family V() f6 (o,ri+r 2 ) defined as 

' /x, 1 if < f < 7i /v* if0<t< Ti 

belongs to C£(0, 7i + r 2 ). 



Mr:=< , v r :=<^ ' (4.14) 



Lemma 4.5 (Compactness for solutions of the continuity equation (I)) Let (fi n ,V n ) be a 
sequence in C£(0, T) such that 

1. for some s e [0, T] sup„ 6N $ (B R ) < +°° VR > 0; 

2. the sequence of maps t h- > v" (S«) w equiintegrable in (0, T), /or every > 0. 
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Then there exists a subsequence (still indexed by n) and a couple (Mr,V f ) 6 C£(0, T) such 
that (recall (14.4b ) 

J U,"^*Mr wea%* i/i Mj c (R fl ') Vr e [0, T] , 
v"^*v wa%*MMioc(M''x(0,7 , );]R rf ). <4 ' L> 

d4. 15b yields in particular 

f <P(Ht, VtW)dt< liminf V? |f ) dr (4.16) 

Jo «T+~ Jo 

/or every sequence of Radon measures y^*y in Mjt„(R rf ), where <P is an integral func- 
tional as in 43.37b . 

Proof Since v" := / r v"df and m" have total variation uniformly bounded on each com- 
pact subset of M. d x [0, T], we can extract a subsequence (still denoted by /x", v") such that 
in Mi oc (R rf ) and v"^*V in Mi 0C (K c ' x [0, T};R d ). The estimate g2) shows that 

su Pj u,"(5s) <+°° Vf e [0,T], /?>0. (4.17) 

The equiintegrability condition on V" shows that V satisfies 

\V\{B R xl) =Jm R (t)dt V/6S(0,T),fi>0, for some Mr e L 1 (0, T) , 

so that by the disintegration theorem we can represent it as V = J Q T V, for a Borel family 
{ v r }te{0,T) st iH satisfying 44.2b . Let us now consider a function £ 6 C c ' (R. d ) and for a given 
interval /= [fo,fi] C [0, T] the time dependent function £(t,x) := X i (t)'V ^ (x) . Since the dis- 
continuity set of £ is concentrated on N = M. d x {to , t\ } and V | (N) = 0, general convergence 
theorems (see e.g. [3 Prop. 5.1.10] yields 



lim / / V£(x)-dV?(x)d/= lim / £-dv"(f,x) 
= f C-dvM= / / VC(jc).dv t (*)dr. 

JH rf x(0,r) J/iR rf 

Applying 44.5b with £(r,x) := f (jc) and fo := s and the estimate 44.17b we thus obtain the 
weak convergence of fj," to a measure ji, e M + (K £/ ) for every t e [0, J]. It is immediate to 
check that the couple (/i, , V, ) belongs to C£ (0, T) . 44.16b follows now by the representation 

J T <P(^,v,\y)dt = <P{n,v\f), n :=y Mr*, 7= 7®^' e M+ C (K J x (0,7)) 
and the lower semicontinuity property stated in Theorem l2.ll □ 
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4.2 Solutions of the continuity equation with finite ^-energy 

For all this section we will assume that <j> : (0, +°°) x K d — > (0, +°°) is an admissible action 
density function as in d3.2b .b.c) for some p € (1,+°°), 76 M 1 J )c (R rf ) is a given reference 
Radon measure, and <P is the corresponding integral functional as in ( 13.37b . We want to 
study the properties of measure valued solutions (jlt,V) of the continuity equation d4.lt with 
finite <f>-energy 

E - / <P(n,,v,\y)di<+°°. (4.19) 
Jo 

We denote by G£* v(0, r) the subset of C£(0, T) whose elements (fi, v) satisfies d4.19b - 

Remark4.6 If (/X/) /e [o,r] ls weakly* continuous in M^ c (R rf ) and ( 14.19b holds, then [i,,V, 
also satisfy ( 14.2b : in fact, the weak* continuity of fi, yields for every R > sup, g j r j fi t (Br) = 
M R < +°°, and the estimate ( 13.48b yields (recall ( 13.43b ) 

Vr<t] r||v,||(5 s )dr<r|r 1 /«Z< 1 /''//(y(5^),M R ) 1 /?<+oo. (4.20) 
Jo 

Recalling that the function h is defined by ( 13.44b . we also introduce the concave function 

™ {s):= IoW)V~« dr > ° ,(0) = ' ffl ' W= h(^' »(*) = +-■ (4.2D 
In the homogeneous case (j)(p,z) = p a ||z||* we have 

CO(s) = dr = S— s l - a li = '- s s Ip. (4.22) 

w Jo q-a 6 

For given nonnegative £ 6 C c 1 (R d ) and ji £ M 1 [ )c (]R fl ') we will use the short notation 

Z := supp(DC) C W, G p (C):=/>d 7 , D(£) := sup ||Df||*. (4.23) 

Theorem 4.7 Le/ £ 6 C ( !(M £i ') be a nonnegative function with Z,G(£),Z)(£) defined as in 
( 14.23b . W Zef ju,v 6 e£^ >y (0,r). Sefft'ng 

E z := { 4>(m,, V, 1 7, Z)dr (4.24) 
Jo 

we have 

\^^ t (i;n\<pD(Q^,Mr,Z) 1/p H{G p (Q,^ P )) l/q - (4-25) 
/« particular, there exists a constant C\ > only depending (in a monotone way) on h,p, T 
such that 

sup ^,(Cn<Q(MC P )+D(QG p (Q l / c 'E l z /l '+D''(C)E z ). (4.26) 

te[0,T] v ' 

Moreover, ifG p (Q > 0, 

\%a>Uk{C p )/G P (C))\ < ^jb *(^ v »\Y,Z) 1/P fara.e.te(0,T). (4.27) 
In particular, in the (a-Q)-homogeneous case, for every < s <t < T we have 

\UC M - llf 11^)1 < t>D(Q IKII^z) /^(M,-,V,| r ,Z) 1 /"dr. (4.28) 
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Proof Setting m, := lk(£ p ), G = G P (Q, D = D{Q we easily have by (3g7) 

jm, = j f J^;''d^=pJj''- l V;-d V ,<pD^ t , Vt \ 7 ,Z) l ^H(G,m t ) 1/q , 

since {^ p ~ l ) q = C p - Sinceff(G,m,) = Gh(m,/G), we get 

h- l /i(m t /G)—m, < P DG 1/q <P(n,, v, \y,Z) 1/p . 

Recalling that ^0)(r) = h _1 / ? (r) we get j4.27> . 

In order to prove ( 14. 26b we set M := sup, g [ r i m r and we choose constants (a, b) e fj, 
integrating d4.25b we get 

1/9 f („r£\ l liv x IP J. (un4\ l liv l IP 



sup |m, - m 1 < pDT 1 '" (aG) L,q E^ p + (bM) l ' q E% p . (4.29) 
re [0,7"] V " ' 



By using the inequality xy < p 1 x p + q 1 y q we obtain 

M < m + /?Z5 (aFG) I/9 /4 /p + -M + p'^ 1 D p {bT) plq E z (4.30) 

which yields fl~26l > with Ci := pmax (l.p^r) 1 /?,^- 1 ^)^). 

Finally, let us assume that <j> satisfies the ( a- 8) -homogeneity condition, so that co(s) = 
\s Q l p as in fi~22l . It follows that 

C0(G-V) = fllCllV)ll£llZ, e ( y,zr (4-3D 

Integrating ( 14.27b we conclude. □ 
We extend the definition of m r (lt) also for negative values of r by setting 



m r (H):=H(Bi)+ \x\ r dpi{x) = 1 V \x\ dit(x) VreR. (4.32) 

JR d \Bi JW 

Notice that rfto(M) = [i(M, d ) and m r (lt) < m r (ll) < + m r (ii) when r > 0. 

Theorem 4.8 Let us assume that m, (y) < +°°/or some r < p and let (li, V) e e£^ ir (0, 7 1 ) 
satisfy ( |4.19t . For every 6" < 1 +r/q, if rfig(tio) < +°° a^o rhg (jU, ) < +°° and r/iere 
exists a constant C2 onfy depending in a monotone way on h,/?, 7", A, |5| swcw fnaf 

m 5 (M/) < C 2 (m 5 (iio) + m r ( r ) l /?£ l /"+ii). (4.33) 

Moreover, ifr > —q and pU)(M. d ) < +°°, then pi, (R d ) is finite and constant for every t 6 [0, T\. 
Proof Let us first set 

K n :=2 nr Y(B 2n+l \B 2 n) (4.34) 

observing that 

+00 

K„ <Y, K J< 2 '' nv(r)> limsupTT,, = 0. (4.35) 

We consider the usual cutoff functions £„ 6 C™(M. d ) as in | |4.1 lal b) and we set 

D„ = D(£„) = sup||D£ n ||» < A2-", G„ = G P (Q < y(5 2 „ + i \5 2 „) = 2-" r K n . (4.36) 
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By ( 14.26b we obtain 

sup ^,{B 2 n) < C l U ) {B 2 ,,+ i )+A2- n(l+r I^Kl lq E x l<i+AP2- n PE)- (4.37) 
te[Q,T\ v ' 

in particular, if r > — q and Ho(R. d ) < +°°, we can derive the uniform upper bound ji, (Mr) < 
C\Ho(M. d ) letting n j +°°. We can then deduce that H t (M. d ) is constant by applying the 
estimate d4.25b . which yields after an integration in time and for every (a,b) 6 fj, 



sup | M ,(Cf)-Mo(Cf)| <pAr 1 /^ 1 ^2-"(a2-"'-^ + foC li Uo(R d )) 



1/9 



«e[o,r 



In order to show ( I4.33I I. we argue as before, by introducing the new family of test functions 
induced by v„(x) := Vq(x/2") e C""*" 



fu„(jc)=l if 2" < \x\ < 2"+\ 

I U„(x) = if |jc| < 2 or |jc| > 2 , 

Observe that 1 < (v„(x)) p < 3 and for some constant Ag > 1 

A s' \ X \ S ^ £2 s "(u„(x)) p <A 5 |;t| 5 VxeR 17 , |jc| >2. (4.39) 



n=l 



As before, setting K' n := K n+ \ +K„-i, we have D(u„) < A2 " and 

G p (u„) < (2-(" +1 ^„ +1 +2-("- 1 ^„_ 1 ) < 2l r l 2-" r ^. (4.40) 
Applying d4.26b we get for every t g [0, T] 

2"V/«) < Cj (2" 6 Mo«) +A2( 5 - 1 - r / f ')"«) 1 /'? {E' n ) l lP +A"2 {S -P>E',}j , 

where 

£„:= / 4>(M,,V,|r,5 2 „ +1 \S 2 „)dr, #;:=£„+!+£„_!. (4.41) 
Since 8 < l + r/q and 5 < p, summing up with respect to n and recalling d4.37b we get 

m 5 ( i u,)<C 2 (m 6 ( i u ) + (m,.(7)) 1 / c 'Zi 1 ^+ J E). □ (4.42) 
In the the 0-homogeneous case we have a more refined estimate: 

Theorem 4.9 Let us assume that is Q -homogeneous for some 6 6 the measure y 

satisfies the r-moment condition rn r (y) < +°°, and let (p.,V) 6 G£^, y(0, T) satisfy d4. 19b . 
For every 8 < 8 := 4/3+ (1 — A)r, ifmg(no) < +°° "is (Mr) is finite and there exists a 
constant C3 > such that 



m 5 ( M ,)<C3(m 6 (/io) + m r (7) 1 - 1/e £ 1 / e ). 



(4.43) 



Moreover, if 8 > (i.e. r > -p/{6 - I)) and fJo(M. d ) < +°° then l± t (WL d ) is finite and constant 
forte[0,T}. 
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Proof We argue as in the proof of Theorem | |4.9I >. keeping the same notation and using the 
crucial estimate d4.28l >. If £„ are the test functions of d4.11al b). 

II & llfT(U) = G «~ l)IP m 2- nr{6 - 1)/p (Kn) (9 - 1)/P , (4.44) 
so that, since 86/ p = 1 + {9 - l)r/p, <4.28t yields 

\{Mtf)) e/P -(Mtf)f P \ <AB2-' ss "IP{K n ) (6 - l)lp E l IP. (4.45) 

Since 8 > 0, passing to the limit as n | °° and recalling d4.35l l. we get Ht(M. d ) = jUo(R rf ). 
Concerning the moment estimate, we replace £„ by u„, defined by in d4.38b . obtaining 

|(M,(<)) e/ ''-(^(<)) e/p |<C3,2- S ^(<)< e - 1 )^(£;0 1/ ^ (4.46) 

and therefore 



Mv<;)<Ci. 2 (M<)+2- dn (Ky ' {K) 1 )■ (4-47) 
Multiplying this inequality by 2" s , summing up w.r.t. «, and recalling d4.39K we obtain 

m 5 ( J ix,)<C 3 (m 6 ( 1 u ) + m,.(r) 1 - 1 / e £ 1 / 9 ). □ (4.48) 



Corollary 4.10 (Compactness for solutions of the continuity equation (II)) Let (fi",v") 

t, + oc( 



be a sequence in C£qj(0, T) and let y"^*y in M ] J )C (K £/ ) such that 



sap fi$(B R ) < +oo Vfl>0, sup / T 4>(/i ( ",v^|y")df<+°o. (4.49) 

neN neN-'O 



TTiew conditions 1. and 2. ofLemma W. 5\ are satisfied and therefore there exists a subsequence 
(still indexed by n) and a couple (jl t ,V t ) G 6£,j>,y(0, T) such that 

tf^*H, weakly* in Mj c (M d ) Vr G [0, T] , 
v"^*V wea%* ;n Mi oc (M d x (0,T) ;WL d ), 

[ T <P(li t ,V t \y)dt <liminf f <P(tf,v';\Y')dt. (4.51) 

Jo «T+°° Jo 

Suppose moreover that Mo(R a ') -> Mo(R rf ) and 

sup„ ^^(y 7 ) < +°° where K — — q or K — 



-p/(9 — 1) in the 8 -homogeneous case, then (along the same subsequence) fJ."(M. d ) 
H,{R d ) for every t G [0,7*]. 



Proof Since Pr := sup,, y"(B R ) < +°° for every R > 0, the estimate d4.33b for 8 = and 
the assumption fl4.49b show that Mr = sup„ 6N te r r j [i"(Br) < +°° for every 7? > 0. We can 
therefore obtain a bound of || V" | by 0.48b . which yields 

\\v 1 }\\{B R )<H{P R ,M R ) llq ^{li,,V,\y) llp , 

so that the maps 1 i— » ||v"||(Sr) are uniformly bounded by a function in L p (0,T). The last 
assertion follows by the fact that 1 1— » /i"(K rf ) is independent of time, thanks to Theorem |4.9l 
(in the (a-0)-homogeneous case) or Theorem l4.8l (for general density functions (j)). □ 
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5 The (0-y)-weighted Wasserstein distance 

As we already mentioned in the Introduction, BENAMOU-BRENIER [7] showed that the 
Wasserstein distance W p ( II. lb can be equivalently characterized by a "dynamic" point of 
view through d 1 . 1 5 b . involving the 1-homogeneous action functional i ll. 13b . The same ap- 
proach can be applied to arbitrary action functionals. 

Definition 5.1 (Weighted Wasserstein distances) Let y e M 1 J )C (M £ ') be a fixed reference 
measure and : (0,+°°) x H d — > [0,+°°) a function satisfying Conditions d3.2b .b,c). The 
(0, y)-Wasserstein (pseudo-) distance between jUo, jUi € M^ c (M d ) is defined as 

W5 )y (Mo,Mi):=inf{/ *(A,V,|y)dr: (/i,v) e 68(0, l;/io -> Mi)}- (5-D 

We denote by M^,y[/io] tne set °f a U the measures ,u e Mj| x .(R' / ) which are at finite W fr - 
distance from /Xo- 

Remark 5.2 Let us recall the notation W 1Kar/ of i ll. 23b in the case p c( (p, w) = p a \v//p a I 7 '. 
When a = we find the dual homogeneous Sobolev (pseudo-)distance i ll. 7b and in the case 
a = 1 and supp(y) = R d we get the usual Wasserstein distance: 

11/4) - Mi =^.o ;r (Mo,Mi), W r p (/io,Mi)=%,i;y(rt) ) Mi)- 
Remark 5.3 Taking into account Lemma |4~31 a linear time rescaling shows that 

W* y fjio,Mr) :=is&{T*- 1 J *(fi,,V t |y)df : (/t,v) e e£(0,r ;j uo - Mr)}- (5.2) 

Theorem 5.4 (Existence of minimizers) Whenever the inflmum in < !5.1b w a finite value 
W < +°°, it is attained by a curve (jU, v) € G£^y(0, 1) smc/i f/iaf 

4>(M,,v,|y) =W forJf l -a.e.t<£ (0,1). (5.3) 

77ie curve (/t;)fe[o,i] associated to a minimum for \5 . 1 b is a constant speed mimimal geodesic 
for "Wa y Ji'wce if satisfies 

Wf, r (H,,jO = l*-'f|'W*,7(W>,Mi) Vs,/e[0,l]. (5.4) 
We /iave a/so f/ie equivalent characterization 

W^, y (CT,Tj) =irf{J (*(M/,v,|y))' /P d ? : (ji.v) 6 C£(0,:T;ct -> r,)}. (5.5) 

Proof When W^^/io, Ml ) < Corollary 14. lOl immediatelv yields the existence of a min- 
imizing curve (,U,v). Just for the proof of l !5.5b . let us denote by Wa ;7 ((7,TJ) the infimum 
of the right-hand side of l !5.5b . Holder inequality immediately shows that W^ t y(a,rj) > 
W^, r (a,r)). In order to prove the opposite inequality, we argue as in [ 3 . Lemma 1.1.4], 
defining for (jU,V) 6 e£(0,r;CT -> 77) 

Se(0 := y' (e + 4>(M r ,v,|y)) 1/P dr, te[0,T}; (5.6) 
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s e is strictly increasing with s' £ > e, s £ (0,r) = (0,S £ ) with S e := s E (T), so that its inverse 
map t £ : [0,S E ] — > [0, T] is well defined and Lipschitz continuous, with 

t £ os £ = (e + $(M,,v,)) ^ a.e.in(0,r). (5.7) 
If /l e = ot £ ,v £ := t' £ Vot E , we know that (£ £ ,v £ ) e C£(0,5' £ ;ct — > r]) so that 

the latter integral being less than S P . Passing to the limit as e J, 0, we get 

W^ r (a,t])<J ^,v,\Y) 1/p dt V( i u,v)ee£(0,r;CT^77), (5.8) 

and therefore W^ r (CT,T)) < W^^ff,?]). If (jU,v) £ e£^,, y (0, 1;jUo — ► /Xi) is a minimizer of 
d57Tb . then ED yields 

W 1 ^ = W^, 7 ( J uo,/Xi)= (^ 1 4>(M,,v,|7)df) 1/ " = ^ 1 4>( i u,,v,|y) 1 /"df, 
so that E3 holds. □ 



5.1 Topological properties 

Theorem 5.5 (Distance and weak convergence) The functional Wi, is a (pseudo)-distance 
on Mjt (R d ) which induces a stronger topology than the weak* one. Bounded sets with re- 
spect to Wi« are weakly* relatively compact. 

Proof It is immediate to check that W^, y is symmetric (since 0(p,— w) = 0(p,w)) and 
W<j> y(cr, 7] ) = <7 = tj. The triangular inequality follows as well from the characteriza- 
tion d5 - 5 b and the gluing Lemma l4~4l 

From ( 14.271 ) (keeping the same notation d4.23l l) and d5 - 5 b we immediately get for every 
Mo, Mi eM+ c (R rf ) and nonnegative ^6C c 1 (R rf ) with ||Clb( 7 ) >0 

\co{M^ p )/G P (Q)-co{^ p )/G P (Q)\ < -^^-W^a.Tj), 

1 1 llslliP(7) 

which shows the last assertion, since CO is strictly increasing and the set 

{p-.CecX® 1 ), C>o, ||CIIlp(7)>0} 

is dense in the space of nonnegative continuous functions with compact support (endowed 
with the uniform topology). □ 

Theorem 5.6 (Lower semicontinuity) The map (mo,Mi) l— * W0,y(Mo,Mi) !S bovver sem;'- 
continuous with respect to weak* convergence in M.^ K {M. d ). More generally, suppose that 
y^*Y in M 1 J )C (R £ '), 0" is monotonically increasing w.r.t. n and pointwise converging to 0, 
and ^^*^q,^^*Hi in Mj c (R d ) as n f +°°. 77iera 

liminfW ri/ .(^,M") > W t , r (Mo,Mi)- (5-9) 
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Proof It is not restrictive to assume that W^n,yn(^ig,/i") < S < +°°, so that we can find a 
sequence (jx n ,v n ) e e£^n t yi(0,l;^ — > jx") such that 

* m ( J U ( ",vf|/') <5 a.e.in(0,l), Vm < » G N, (5.10) 

where <f>"' denotes the integral functional associated to (j) m . We can apply Theorem |4. lOl and 
we can extract a suitable subsequence (still denoted bu fi",V n ) and a limit curve (jl,v) 6 
e£0, r (O, l;/Xo — » Ml ) sucn mat < !4.50b holds. We eventually have 

W'; m , 7 (Mo,Mi) < j 4> m (ti„V t \y)dt<S. (5.11) 
Passing to the limit w.r.t. m | +°° we conclude. □ 

Theorem 5.7 (Completeness) For every a eM+ c {R d ) the space M^.y[ff] endowed with 
the distance Wi« u complete. 

Proof Let (p.„)neN be a Cauchy sequence in M^.y[(r] w.r.t. the distance W^.y; in particular, 
(fi„) is bounded so that we can extract a suitable convergence subsequence \i n . weakly* 
converging to fj. x in MjJ XJ (R d ). Thanks to Theorem 15.61 we easily get W ( (,,y( J u m , J u 00 ) < 
liminft^ 00 "W^.y( J u m , J u„ (: ), and therefore, taking into account the Cauchy condition, 
limsup J „^ 00 W^, T ( J u m , J u 00 ) < limsup n m ^ oo W 0j y( J u m ,jLi„) = 0sothat,u„ converges to □ 

Let us now consider the case of measures with finite mass (just to fix the constant, probability 
measures in 5 > (R rf )). We introduce the parameter 

1 if <j) is (a-0)-homogeneous, 



0-1 1-a " (5.12) 

FT 



-2— — q otherwise. 



Theorem 5.8 (Distance and total mass) Let us assume that m- K (y) < +°° and let us sup- 
pose that o~ € 7{R d ). Then M* y[o\ C 7{M. d ), the weighted Wasserstein distance W^, v is 
stronger than the narrow convergence in T(M. d ), and 7(M. d ) endowed with the (pseudo-) 
distance VJ^.y is a complete (pseudo- )metric space. 

Proof If 7] 6 M^.yfcr] then Theorem 14.91 (in the 0-homogeneous case) or 14.81 (in the gen- 
eral case) yields tj (R d ) = a(M. d ) = 1, so that M^ r [a] c 9(R d ). Since the narrow topology 
coincide with the weak* one in 3 > (R £ '), Theorem |5.5| proves the second statement. The com- 
pleteness of !P(R rf ) with respect to the (pseudo) distance W<j,,y follows by Theorem l5.7l □ 

We can also prove some useful moment estimates. 

Theorem 5.9 (Moment estimates) Let us assume that rh, (y) < +°°for some r e K and let 
us set 

g ,_ jip+C 1 - i) r = f ( l + r / K ) if $ is -homogeneous, 
\l + r/q<p otherwise. 

If m§(<y) < +°°for some 8 < 8, and 7] £ M* y[ff], then rr\$(Tj) is finite and there exists a 
constant C only depending on (j),8 such that 

m 5 (7])<c(m 5 (0-) + m,(y)+W'; /ff.Tj)) 

; , (5.i4) 

m s (7 7 )<C(m 5 (cT) + m r (r) 1 " 1/e ^ / a e (a,r ] ) 



Moreover, when 8 > 1, the topology induced by "W^y in Mi v[cr] is stronger than the one 
induced by the Wasserstein distance Wg. 



2S 



Proof Let us first consider the general case: applying d4-.33b we easily obtain d5. 14b . In 
order to prove the assertion about the convergence of the moments induced by "W^y (which 
is equivalent to the convergence in Wg when 8 > 1), a simple modification of ( 14.42ft yields 

/ |*| 5 d77<C 3 (7 W 5 da + m r (7) 1 /*W f7 (a,77)+< y (a,T])), (5.15) 

J\x\>2" vj|.t|>2"- 1 v '" ' 

which shows that every sequence r\ n converging to <7 has 5-moments equi-integrable and 
therefore it is relatively compact with respect to the 5-Wasserstein distance when 8 > 1. 
The 0-homogeneous case follows by the same argument and Theorem l4.9l □ 

Remark 5.10 There are interesting particular cases covered by the previous result: 

1. When 7(K £/ ) < +°° then W^.y is always stronger than the 1-Wasserstein distance W\ ; in 
the ^-homogeneous case, W p , a ^ also controls the W p jq distance. 

2. When m p (y) < +°°, then W<j,,y is always stronger than W p . 

3. When (j) is ^-homogeneous with 9 > 1 and y is a probability measure with finite mo- 
ments of arbitrary orders (this is the case of a log-concave probability measure), then all 
the measures a € M^.yfy] have finite moments of arbitrary orders and the convergence 
with respect to W^.y yields the convergence in CPg(R << ) for every 8 > 0. 

5.2 Geometric properties 

Theorem 5.11 (Convexity of the distance and uniqueness of geodesies) (•,•) is con- 
vex, i.e. for every \l\ <E Mj+ C (R d ), i, j = 0,l, and T £ [0, 1], if ^ = ( 1 - t) M° + Tflf , 

w; r ( Mo T X) < (i-T)w; r ( Mo °,/i?)+TW'; r ( M( j, Ml 1 ). (5.i6> 

If (j) is strictly convex and § has a sublinear growth w.r.t. p (i.e. = 0), then for ev- 
ery Ho, Hi € M 1 J )c (R fl ') with W^.y(/io,A I i) < +°° there exists a unique mimimizer (/X,v) € 
e£ f y(0,l)o/(EB- 

Proof Let (M^V- 7 ) £ C£^. 7 (0, l;Mo — ► Mi) be two minimizers of d5.1l l. 7 = 0,1. For re [0, 1] 
we set/i, T := (1 -t)m, + tm/, vj := (1 - t)v ( ° + xvj. Since (m t ,v t ) £ C£(0, l;/x T -> M?). 
the convexity of <j> yields 

wy^oX) < [*Wtf\Y)*<f ((l-T)*(tt°,v?|y) + T*OV,v/|y))df 

= (i-T)w'; 7 (< / i 1 °)+Tw; 7 ( / i 1 , Ml 1 ). 

Let us now suppose that <j> is strictly convex and sublinear. Setting, as usual, fif = p f T y+ 
(jU^) 1 , vj = w f T y, we have for a.e. t 6 (0, 1) 

4>(M, t ,v, t | 7 )<(1-t) f *(p ( °,w?)dy+T / ^(p/.w^dy (5.17) 

and the inequality is strict unless p r ° = p, 1 and = w, 1 for y-a.e. x g M £/ . If Mo = Mo anc ^ 
Ml = Mi' two minimizers should satisfy 

p°(z) =p/(x), w?(x) =w/(x) y-a.e., V? = V, 1 for J^'-a.e. f e (0, 1). 
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Since (m',V) are solutions of the continuity equation, taking the difference we obtain 

d, ((M,°) X - (M/ ) X ) = d, (M,° - M, 1 ) = -V ■ (v? - v, 1 ) = in R d x (0, 1). 

The difference (Mr°) X ~~ (m/) X ls tnen independent of time and vanishes at t = 0, so that 
M,° = jll/ for every t £ [0, 1]. □ 

Theorem 5.12 (Subadditivity) For every m/ e M^ oc (R d ), i,j = 0, 1, we have 

W f r (/x ° + Mo , Mi° + Mi 1 ) < W^, r (Mo°, + W^, r (^ , Mi)- (5.18) 

In particular 

W^ r ( Mo + a, i Ui+cr)<W^ r ( J uo, J Ui) V<7GM+ c (M d ). (5.19) 
Proof Let (m^V 7 ') 6 C£,/>, 7 (0, l;Mo -> Mi) be 

as in the proof of the previous Theorem. Since 
( J u°+M 1 ,v°+v 1 ) e e£(0,l;Mo+Mo ^M°+Mi 1 )^ we get 

W^, 7 (Mo° + Mo' , M? + Mi 1 ) < jT 1 (<*>(M,° + M, 1 , V? + v, 1 |y)) V " dr 

< jT' [(€>(^ + ^^\y)) 1/P +(<P(^ + ^,vj\y)) l/P ]dt 

< jT 1 [(4>(/x,°,v; ) |r)) 1/P + (^(M/^/lr))' 7 "] df = W f7 (M^M?)+W^ 7 C<Mi)- □ 



Proposition 5.13 (Rescaling) For every Mo, Mi £ Mio C (KL d ) and A > we have 

w; a7 (Amo,Ami) = AW; 7 (mo,Mi), (5.20) 

W^ r (AMo,AMi)<A^ r (Mo,Mi) «TA > 1 
^ 7 (AMo,AMi)<Aw; y (Mi,Mi) if A < 1. 

Proo/ J5.20b follows from the corresponding property <P(Xli, Av|Ay) = X<P([i,v\y). Anal- 
ogously, the monotonicity and homogeneity properties of <j> yield 

0(Ap,Aw) <0(P,Aw) = A p 0(p,w) if A > 1; 

the convexity of <j> and the fact that <j> (0, 0) = yield 

0(Ap,Aw) < A0(p,w) if A < 1 . 

( 15.2 1 b follows immediately by the previous inequalities. □ 

Proposition 5.14 (Monotonicity) Ify\>yi and 0i < fa then for every Mo, Mi e ^ioc(^) 
we have 

W^. 7i (mo,Mi) < W fcjR (Mo,Mi)- (5-22) 



Theorem 5.15 (Convolution) Le? k 6 C~(R ) fee a nonnegative convolution kernel with 
(x)Ax= 1 and letk e (x) := e~ d k(x/e). For every Mo, Ml £ •M-iocC® ) we ^ ave 

W^, rrfE (Mo*^,Mi *^e) < W^,. 7 (Mo,Mi) Ve>0; (5.23) 
limW^, rrfE (Mo*fc£,Mi *^e) = W^,. r (M),Mi). (5-24) 
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Proof Let (jll,v) 6 C£,/ )) y(0, l;^o — ► Ml) be an optimal connecting curve as in Theorem l5.1 II 
and let us set \if = ^*k e ,vf := v r *k £ . Since (pi £ ,V £ ) e G£(0, 1;/Xq 45.23b then 

follows by i2. 12b whereas ( 15. 24b is a consequence of Theorem 1 5. 6 1 □ 

Remark 5.16 (Smooth approximations) For a given curve (ju, v) € G£f y(0, 1) the convolu- 
tion technique of the previous Theorem exhibits an approximations (fi E ,v £ ) in ye (0, 1), 
~f := y*k £ , which enjoy some useful properties: 

1. p: £ = p £ .if d , V £ = w £ ^f d with p e ,w E e C°(R d ); if Mo(K fl ') < +°° and m_ K (y) < +°° 
(recall Theorem |5. 8b , then p £ is also uniformly bounded. 

2. If supp(£) C ~B\ then p e ,w e are supported in G £ := {x eR d : dist(jc,G) < e}, G = 
supp(y). 

3. p e ,w £ are classical solution of the continuity equation 

d,p £ + V-w £ = inK d x(0,l). 

4. If (/J.,V) is also a geodesic, / (j)(p £ ,wf)df < <&(p t ,V,\y) = v (Mo,Mi) for ever y 

re [0,1]. 



5.3 Absolutely continuous curves and geodesies 

We now study absolutely continuous curves with respect to "W^.y and their length. Let us 
first recall (see e.g. [3 , Chap. 1]) that a curve (h>/j,£ Mi oc (K £/ ), t e [0, J], is absolutely 
continuous w.r.t. W^y if there exists a function m £ ^'(0, T) such that 

Wf, r (|^,rtb)< f'm^dt V0<t <h<T. (5.25) 



The curve has finite p-energy if moreover m € L p (0,T). The metric derivative |/x'| of an 
absolutely continuous curve is defined as 

|K|:=li m %^, (5.26) 

h^Q \h\ 

and it is possible to prove that |/x/| exists and satisfies |/x/| < m(f) for Jzf ! -a.e. ? 6 (0, J). The 
length of jU is then defined as the integral of \fi'\ in the interval (0, T). 

Theorem 5.17 (Absolutely continuous curves and their metric velocity) A curve 1 1— > \i t , 
t 6 [0, T], is absolutely continuous with respect to W^.y if and only if there exists a Borel 
family of measures {Vt)te{o.T) m Mi oc (]R £ ';R rf ) such that (jU,V) 6 e£^, i7 (0,r) aw/ 

jf (4>(M„V,|r)J df<+-. (5.27) 

In this case we have 

\l4\ p < <P{n,,v,\y) for£' l -a.e.te(Q,T), (5.28) 
and there exists a unique Borel family V t such that 

= *(M,,V»|y) for^-a.e. t e (0,T). (5.29) 
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Proof One implication is trivial: if (jU,v) € C£^, r (0, T) and ([5727} holds, then d5.5l l yields 

W^ fl ,M f0 ) < jf* (^(Mf.^lr)) 1 V (5.30) 

so that fi is absolutely continuous and d5 .28b holds. 

Conversely, let us assume that pt is an absolutely continuous curve with length L. A 
standard reparametrization results (3] Lemma 1 . 1 .4] shows that it is not restrictive to assume 
that fj. is a Lipschitz map. We fix an integer N > 0, a step size T := 2~ N T, and a family of 
geodesies ,V k>N ) <E e£A j7 ((fc- l^&Tj/J^w — > jX kt ),k= 1, - ■ ,2^, such that 

T*(M^,v*^|y) = t 1 -'^ J^dt,/^) < r dr. (5.31) 

J(k-\)x 

Let (n N ,V N ) e C£ A, y (0, T) be the curve obtained by gluing together all the geodesies (fi k ' N , V k 
Applying Corollary 14.101 we can find a subsequence (n Nh ,V Nh ) and a couple (jlt,v) £ 
e£^, r (0,r) such that fl^^fl, for every f e [0,7] and v^^V in Mi oc (R'' x (0,7 , );R £/ ). 
It is immediate to check that fa = fa for every t 6 [0, T] and 

[ T \pi' t \ p dt > liminf C ' ®{\£ h ,vf*|y)df > f <&(/x,,V,|y)d* > f |A,'| p d/, 
Jo &t+°° Jo Jo Jo 

which concludes the proof. □ 

Corollary 5.18 (Geodesies) For every fj. e M 1 J )c (R' i ) f/ie space Ma y [jll] is a geodesic space, 
i.e. every couple Ho, fa 6M^.y[jU] can be connected by a (minimal, constant speed) geodesic 
t e [0, l] jU, e MA. y [/z] jwc/i f/za? 

W^, r (M,,/i f ) = |*-*|'W^ r (/io,Mi) V*,*e[0,l]. (5.32) 

A// f/;e ( minimal, constant speed) geodesies satisfies the continuity equation l l4.lt for a Borel 
family of vector valued measures {v t ) te (p n such that 

*(Mf,v,|r)=W5 )y (Mo,Mi) /We. re (0,1). (5.33) 
//" is strictly convex and sublinear, geodesies are unique. 

Remark 5.19 (A formal differential characterization of geodesies) Arguing as in (27] Chap. 
3], it would not be difficult to show that a geodesic fa = P/-S? rf with respect to W 2 a .cgd 
should satisfy the system of nonlinear PDE's in R d x (0, 1) 

a „ ,,_ ,2 
dtY+^P |Vy| =0= 

for some potential t//\ Unlike the Wasserstein case, however, the two equations are coupled, 
and it is not possible to solve the second Hamilton-Jacobi equation in \jf independently of 
the first scalar cosnervation law. In the present paper, we do not explore this direction. 
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We can give a more precise description of the vector measure V satisfying the optimality 
condition l |5.29l l. For every measure /i 6 M 1 [ )c (K rf ) we set 

Tan 0j7 (/i) := {v e M loc (R c ';R c ') : #(/i,v|y) < +°°, 

, , 1 (534) 

*(M,v|y) < 4>(M,v + Tj|y) Vtj e Mioc^iM 17 ) : V-7] =0|. 

Observe that for every v € Mi oc (K rf ;M d ) such that #(/x,v|y) < +°° there exists a unique 
V :=IT(V) e Tan^,, r ( J u) such that V- (v-v) =0. In fact, the set AT(v) := {v' e Mi oc (IR fl ';IR fl ') : 
V(v' — v) =0} is weakly* closed and, by the estimate d3.47l > the sublevels of the functional 
v' i— » <P(fi,v'\y) are weakly* relatively compact. Therefore, a minimizer V exists and it is 
also unique, being 4>( J u,-|y) strictly convex. 

Corollary 5.20 Let {pL.v) e C£,j,. r (0, T) so that [1 is absolutely continuous w.r.t. W^.y. The 
vector measure V satisfies the optimality condition 1 15.291 ) if and only ifV, 6 Tan^ ,y(/i,) for 
& l -a.e.t&(Q,T). 

Let us consider the particular case of Example 13. 5l in the case of a differentiable norm ■ 
with associated duality map j\ = D|| ■ ||. We denote by j p (v/) = ||w|| p_2 y'i (w) the p-duality 
map, i.e. the differential of ^ || • \\ p and we suppose that the concave function h : [0, +°°) — > 
[0, +°°) satisfies 

\imh(r) = lim r~ l h(r) = 0. (5.35) 

For every nonnegative Radon measure p. 6 M 1 [ )C (]R"') whose support is a subset of supp(y), 
we define the Radon measure h(p.\y) by 

h(p.\y) := h(p) ■ y where p:=^. (5.36) 

Observe that h(fi\y) <C y even if fj. is singular w.r.t. y. 

Theorem 5.21 Let fi e M 1 + c (M d ) and as in d3.23b with h satisfying d5.35l l. A vector mea- 
sure V satisfies 0(jU,v|y) < +°° iffv = \h{p.\y)for some vector field v e L P h , , >(R d ;R d ). 
Moreover, V £ Tan^,.y(/i) if and only if the vector field v satisfies 



rll, 



^(▼)e{vf :f ecrOR'')} *^ ' (5.37) 

Proof Being /; sublinear, the functional 4" admits the representation 

#(M,v|y)= [ J(p,w)dy= I : h( P )\\w/h(p)\\Pdy= f ||vf dA(/i|y), (5.38) 



where /I = py+jU x and v = wy = h(p)\y. The condition V = v/i(/x|y) e Taity, y (/x) is then 
equivalent to 



IvrdM/ily) < / \\\ + z\\ p dh{n\y) Vz e Z/(/*(/i|y)) : V- (zh(u\y)) =0. 
Thanks to the convexity of 1 1 ■ 1 1'\ the previous condition is equivalent to 

J^j p (v)-zdh(n\y)=0 Vze^ (M|/) (M d ;M rf ): J^z-V;dh(p\y) =0, (5.39) 
i.e. /p(v) belongs to the closure of {V£ : £ € C?(M. d )} in L^ |r) (R rf ;K rf ). □ 
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5.4 Comparison with Wasserstein and W Lp distances. 

Theorem 5.22 //y(K rf ) < +°° then for every Ho, me M^ oc (R d ) and a < 1 we have 

Wp/gi^m) < W„/e,uM),Hi) < r(M. d ) 1/K W p , a , 7 (^,fii), (5.40) 
where, as usual, 6 = (1 — u)p + a. 

Proof Let (/X,v) e e£^ pQ ,. 7 (0, 1*,jUq — > /Xi) be an optimal curve, so that 

P-.-fn„ n,\ — [ ch ..(„. \t.\At— [ [ (r,.\ e -P\™.\P 



K«:vU*>,Hi)= **«(Mf,v,)df= / / (p,) B -'|w,|'dydf, (5.41) 

JO JQ JR d 

where p., := p f y+ fJ-i 1 , V, = w,y <S y. Holder inequality yields 



w P /e,i 



. Y 0k>,in)< t I (ft) 1 - p/e |w,| l ' /fl dy*<y(R ,, ) 1 - 1 / fl ^ r 0io,/ii). □ 



Theorem 5.23 Let us suppose that m_jt(y) < +°°, K = p/(d — 1) = #/(l — a), Zei jUo,jUi 6 
3 3 (R c! ), and let K* = k/(k — \ ) be the Holder's conjugate exponent of K. Then 

llMo-Mill^-i.** =W K .,o ;y (/Jo,Mi) < Wp, a:r (Mo,Mi)- (5-42) 

Proof We keep the same notation of the previous Theorem, setting 

1 _ P~ B 



x = p/r:=l+p-e, t*:=— - = 1 + -, x={x*) 



T-l p-6 v 7 l+p-d 

Observing that pi t g J > (R fl ') thanks to Theorem |4.9l we obtain 

W^ r (Mo,Mi)< t />,| r dydf = f f {p t ) x {p t y x \yv t \ r dydt 
' " Jo Jn. d Jo Jr j 

< [ (^p- tT |w,|"dy) 1/T df = (fc j 1tl P - p \«t\ p *rta) l,% = w;, a;r ( Mo ^i). □ 

Theorem 5.24 (Comparison with W p ) Assume that ye MjL(R d ) is a bounded pertur- 
bation of a log-concave measure (e.g. y = f e - v 5£ d , where V is a convex function and f 
nonnegative and bounded). If Pi = SjJ 6 !P(K ') with Si 6 L°°(y) and m / ,( J U,) < L < +°° then 
Wty y(p4), jUi) < +°° and there exists a constant C only depending on L, (j), and y such that 

Wt, r OiQ,Hi)<CW p Qio,Hi). (5.43) 

Proof It is not restrictive to assume that y is log-concave. We can then consider the optimal 
plan E 6 l M + (M. d x R d ) induced by the p-Wasserstein distance dl.ll l between po and pi and 
the interpolant p., defined as 

p,(A)=E({(xo,x l )eR d xR d : (l-t)x +txi e A}) VAeB(R fl '). (5.44) 

It is possible to prove (see e.g. J3] Theorems 7.2.2, 8.3.1, 9.4.12]) that p t is the geodesic 
interpolant between Hq and p\ , it satisfies the continuity equation 

d,p,+V-V, =0 inR d x(0, 1) 



34 



with respect to a vector valued measure V, = \ t ft -C ft where the vector field \ t satisfies 

/ <P p {fi,,v,)di = S [ \y t (x)\ p dft(x)dt = W p (ft),ft), 
Jo Jo JW 

and finally ft = s t y with \\s t \\l°°m <L := max (H^olli^fy); Iki IIl°°(7))- Observe that, being 
s,(x) < L for y-a.e. xeM. d and 0(0,0) = 0, Theorem [TT1 yields 

(j)(s t ,s,\ t ) < ~(j>(L,Lvt) < C L s,\\,\ p y-a.e., 

so that 



5>(ju f ,v,)= / <l>{s„s t Y t )dy(x)<C L \y t \ p s t dy=C L \\,\ p dft, 

Jw d J«. d jR d 

and therefore 

W p Jno,m)< f 1 <P(Ht,v t )dt<C L [ l / \v t \ p dftdt = C L W p {ft>,ft). □ 
v,/ Jo Jo Jm< 

Corollary 5.25 If ft = SjJf d £ 7(M d ) have L°" -densities with compact support (or, more 
generally, finite p -momentum), then ^d(ft), ft) < +°°. 

Theorem 5.26 If fa = s,y with s, > L > y-a.e. in M. d , then there exists a constant C de- 
pending on L and (j) such that 

W f7 (Mo,Mi)<C L ||^-Mi||^i, P - (5.45) 



Proof Let us first observe that if \\ft) — ft ||^-i,j> < +°° then there exists w € L p 



Plod.-m,d\ 



such that 



-V-v = ft-ft h V:=wy, f \yv\ p dy=\\ft-ft\\ p (5-46) 

Jk<i w 7 

In fact, in the Banach space X := Ly(B. d ;M. d ) we can consider the linear space Y := {D£ : 
£ 6 C\ (R d )} and the linear functional 

(t,y):= [ Cd(Mi-Mo) ify = D£ for some £ € C> (R rf ). 

jR d 

lis well defined and satisfies |(£,y)| < \\ft) — P-i W^-^p l|y|lz«(R<'-R<') for every y e 7. Hahn- 
Banach Theorem and Riesz representation Theorem yield the existence of w 6 L p (M. d ;M. d ) 
such that (^,y) = J rJ w y dy, which yields ( 15.46b . Setting ft = (1 - f)/x + fjUi, it is then 
immediate to check that (ft,v) e C£(0, l;/Xo — ► jtii); we can then compute 

W5 r (Mo,Mi)< / / *((l-0*o + ^i,w)dyd/< / A(L,w)dy< C L / |w|"dy, 

where we used the fact that (1 — t)so + ts\ > L y-almost everywhere and the map p i— > 
(j) (p , w) is nonincreasing. □ 
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5.5 The case y = Z£ d and the Heat equation as gradient flow 

One of the most interesting cases corresponds to the choice 

y:=J? d , h a (p):=p a ,0<a<l, $ p . a (p,w):=p a \w/p a \". (5.47) 
In this case the expression of the weighted Wasserstein distance becomes 

H, = p,& d + n t ± , M| ;=0 = W), M|, =1 =J"i}- 
The metric W a .^>d restricted to T(K' i ) is complete if <i < K = = jr^- 

Remark 5.27 CP(M. d ) is not complete w.r.t. W p a -jgd if d > k) The above condition is almost 
sharp; here is a simple counterexample in the case d > K. We consider an initial probability 
measure with compact support Ho = Po^£ d , Po £ L°°(M. d ), and, for t > 0, the family 

H Pt(x):=e- dt p (e-'x), V, :=x\i t =xp t (x)J? d . (5.48) 

It is easy to check that (jU,v) G C£(0,+oo), jU,(R d ) = 1. Evaluating the functional <P, := 
<&p ia (jll„V t |jz? d ) we get 

«k= I A pf- p \p,x\ p Ax= [ e- de 'p^(e- d 'x)\x\"dx 

= e dt-det+pt I e -dt p 8^-d, x ^ e -d tx y, Ax=e (d(l-6)+ P ) t I p{](y)\y\F dy 

Jm.'i jR d 

so that 

^ =ce (i-d/^ ) /" + °° 4)/ i/P df = c _!^ <+00 if d>K . 
Jo d—K 

\f d > K we obtain a curve * 1— > jUf G !P(R d ) of finite length w.r.t. W^, ^.^j (in particular 
(Mn)«6N is a Cauchy sequence) such that lim fit = in the weak* topology. 

In the remaining part of this section, we want to study the properties of W p a .^j with respect 
to the heat flow. We thus introduce 

and for every € M 1 J )c (R fl ') with mg(ju) < +°° for some 5 < 0, we set 

S ; [ J u]=M*ft = "^, «rW=S<MW:= [ g,(x-y)dn(y). (5.49) 

It is well known that u G C°°(K d x (0, +00)) and 

(9,m-4m = inR d x (0,+oo), S t [ju]^*jU as f J. 0. (5.50) 
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Theorem 5.28 (Contraction property) Let ii ,^ 1 € Mj c (IR d ) with m 5 (ju') < +°o 

^/(li 1 , li 2 ) < +°°. T/'jLi/ := 8? are the corresponding solutions of the heat flow, then 

W^M/^fCW^GuV) Vr>0. (5.51) 

Proof It sufficient to approximate the Gaussian kernel g by a family of C°° kernels V with 
compact support and then apply Theorem 15. 151 observing that k" * 5£ d = J?f d . □ 

We consider now the particular case of the W 2 a .cgd weighted distance with a > 1 — 2/d. 
Let us first introduce the convex density function (recall i ll. 9b ) 

^ p):= (2- g) 1 (i- g ) p2 " a - suchthat < {p) =wr p ~ a ' (5 - 52) 

and the corresponding entropy functional 

W a (n)=W a (^\Jf d )-- [ w a (p)dx, ifit = pjSf d <J*? d . (5.53) 



We also introduce the set Q := {/x e 7{R d ) : f (it) < 

Theorem 5.29 If H e 7{R d ) then it, = §,[lt] = 6 Q for every t > 0, the map t 

¥a(/X ( ) is nonincreasing, and it satisfies the energy identity 



f a (M,) + J &i ta (u r ,Vu r )dr = W a {n s ) V0 < s < t < +00; (5.54) 

when jj, 6 Q r/?e« f/ie previous identity holds even for s = 0. Moreover, jl t satisfies the Evo- 
lution Variational Inequality 



2df 



wZ a .jzM,o)+%(Ht)<1'a(o) Vf>0, VaeQ. (5.55) 



Proof Since V*a( M ) = a , a direct computation shows 

d f \jf a (ui) dx = - ^- f Vu, -V\/ a (u,)dx = [ \Vu,\ 2 u; a Ax = 4> 2 .a{u,,Vu t ) 



At Jwi i t 

Concerning d5.55b . we use the technique introduced by 1121 §2]: we consider a geodesic 
(Os,V s )j 6 [q ) i] e C£(0, 1;ct — > it), which satisfies a s (R d ) = 1 by Theorem[4~9l We set 

ct s £ , = ul,J? d := $ £+sl [o s ] , v £ s ,, = yvl,J? d := S e+St [vj , wf, := w s £ , - fV^ . 
It is not difficult to check that 

d s u E s t + V ■ , =0 in R ' x (0, 1) , (5.56) 

so that 

W L;^> (Ve+< , a) < f o A £ SJ As, A% := «,) ~ a | w* ( | 2 dx = *2. a (of, , v*, | J^) . 
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A?, = / (rf,)- 1 -2tVu £ sr w £ st + \w £ st \ 2 -t 2 \Vu £ st \ 2 )dx 



We thus evaluate 

V =/,«*)""(■ 

< -2td s W{af tt ) + ^ a {a s ,v s \^ d ), (5.57) 
where we used the facts 

/ (Mf,) _0! |wf,| 2 dx = 4>2,a(CT s *g £+J/ ,V J *g £+s ,|if c ') < ^ B (o- s ,V s |^f d ) 

JR'' ' 

thanks to the convolution contraction property of Theorem 12.31 Integrating l !5.57b with 
respect to i from to 1 and recalling that (a s , V s ) se p ^ is a minimal geodesic and that 
af , = He+t and CTq , = a, we get 

l A £ ,ds + 2tW a (n £+t ) < 2ff a (a)+W 2 2 a; ^,( J u,CT). (5.58) 



We deduce that 

k^ta ,o)+ tV(Lk+, )<tV(o) + ^wi >a (ii,o). (5.59) 

Passing to the limit as e J, and then as t J. after dividing the inequality by t we get d5.55b 
at t = 0. Recalling the semigroup property of the heat equation, we obtain l !5.55b for every 
time t > 0. □ 

( 15.55b is the metric formulation of the gradient flow of the (geodesically convex) functional 
fa in the metric space (Q, W 2 a ^d), see (3] Chap. 4]. Applying 1121 Theorem 3.2] we even- 
tually obtain: 

Corollary 5.30 (Geodesic convexity of f a ) Let a > 1 - 2/d, M = Pi% d £ V(R d ) with 
W 2 a .^d(fio,Hi) < +°° and / K< / pf~ a dx < +°°, and let fj.,, t £ [0, 1], be the minimal speed 
geodesic connecting Mo to Mi w.r.t. W 2 a -^J- Then for every t 6 [0, 1] jl t = Pt^tf d <C 5£ d 

f p 2 - a Ax< (1-f) / p 2 - a dx + t f p 2 - a dx. (5.60) 

JW 1 JR d JR d 
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